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ABSTRACT 


When investigating the detection performance of a passive 
homing torpedo used against shallow draft surface ships, 
certain environmental factors such as the rough sea surface 
and the bubble-dominated inhomogeneous layer near the sea 
surface have to be considered. This thesis attempts to gain 
some insight into the behavior of a homing torpedo system 
during its critical attack phase, as well as getting some 
indications of the relative importance of the scattering 
mechanisms and the induced tactical limitations. An ideal- 
ized propagation model was used as reference of comparison. 
For a given sea state and target speed the results stress 
the importance of low operating frequency as well as a high 
maximum turn rate. They also point to the importance of 
having a search depth below the bubble-dominated subsurface 
layer, and a variable speed capability during the torpedo's 


attack phase. 
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I. INTRODUCTION 


The following analysis examines several factors that 
limit the detection performance of a passive homing torpedo 
with the mission objective of countering shallow-draft 
targets in Norwegian coastal waters. Generally, these 
factors can be divided into three main groups: 

-Environmental factors in the ocean 

-Electrical mechanical and hydrodynamical factors in 

the torpedo system 

-The users tactical situation 
The factors that are generated in the ocean itself are the 
principal subject of this analysis. As our interest is 
confined to the layer immediately below the surface of the 
ocean, later called the subsurface ocean layer, the main 
Factors affecting the sound propagation are: 

-Scattering and absorption due to the bubble-dominated 

inhomogeneous subsurface ocean layer. 

-Scattering from the rough sea surface 
The concentration of air bubbles and the roughness of the 
sea surface are determined by the windspeed. The effective- 
ness of both these scattering mechanisms depends on the 
frequency of the incident wave and the geometry of the 


source and receiver. The following analysis is limited to 
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high frequencies in the region of 30-60 kHz which are 
characteristic of existing torpedo systems. At high fre- 
quency and low grazing angles for the incident and 
received signals, the phenomena of "shadowing" of the 
surface by other parts of the boundary occurs. Under 
these conditions, the effect from the inhomogeneous sub- 
surface ocean layer becomes increasingly important. Ob- 
viously both the above mentioned scattering mechanisms are 
present simultaneously. Often these two effects cannot be 
resolved either theoretically or experimentally, as any 
Signal with a finite duration will be scattered from the 
space near the surface simultaneously with that from the 
sea surface itself. 

In order to adequately describe the scattering mech- 
anisms, this analysis starts with a presentation of the 
oceanographic background material for: 

-Typical windspeed and wave height 

-Typical sound speed profiles 

-Density and distribution of air bubbles in the sub- 

surface layer 

-Statistical description of the sea surface 
The analysis proceeds by separately estimating the effect 
(eye 

-Scattering from a randomly rough surface 

-Scattering and absorption caused by an inhomogeneous 


subsurface layer, 
13 
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and comparing their relative importance. The method 
employed for these estimations is an approximation that 

1s a combination of both ray and wave theories. Ray 
methods are used to follow the acoustic signal from its 
Source to the vicinity of the scatterer. Wave theory is 
used to calculate the actual scattering process. Finally, 
may theory is again used to follow the scattered signal 
to the receiver. 

An idealized propagation model consisting of an iso- 
tropic stratified medium will be used as reference of 
comparison. This model is founded on: 

-A noise source from a cavitating propeller blade. 

~The operational characteristics for a square law 
detector (ROC-curves). 

-A transmission loss model based on geometrical 
Spreading and absorption losses in homogeneous sea 
water. 

In this analysis, the passive sonar equation is used to 
predict the performance of the homing system. The detec- 
tion range encountering the two scattering effects will 
be obtained from the sonar equation and compared to the 
detection range based on the reference model. Thus, the 
difference in ranges at which the homing device just 
acquires the target with and without scattering 1s a 


measure of efficiency. 
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Ii. SCENARIO AND TORPEDO RUN GEOMETRY 


In the Norwegian coastal waters, the primary mission of 
a torpedo system is to counter an amphibious force consist- 
ing of escorts, supply ships, and shallow-draft landing 
crafts. Typical characteristics of these three ship 
types are as follows: 
Supply ships: 
-Displacement 5000 tons 
-Length 100 m. 
-Draft 6 m. 
=Spe@ed 15-20 kts. 
Escort ships: 
-Displacement 2000-3000 tons 
-Length 85 m. 
-Draft 3 m. 
-Speed 35 kts. 
Laneang craft: 
-Displacement 1000 tons 
-Length 80 m. 
=—DPratt 2 mM. 
-Speed 18 kts. 
In order to simplify this analysis, moderate sea states (SS 3) 


are assumed. Since "moderate" wave heights of two meters 


LS 





are appreciable when compared with the two-meters draft 

of the landing craft, the possibility of an acoustic torpedo 
impacting the target at a depth of two meters is very remote 
without the use of an influence exploder. 

Two relevant search and attack schemes will be con- 
Sidered. These are illustrated in Fig. 1 together with the 
operation of the influence exploder. From target valida- 
tion to completion of terminal attack, the torpedo contin- 
uously tracks in the azimuth plane. In Case A, ascent is 
inhibited after enable. For Case B, ascent is inhibited 
after the torpedo reaches terminal attack depth. 

An assumed attack depth of six meters is consistent 
with the activation range of influence exploders and is 
deep enough to preclude wave or "free surface" induced 
disturbances of the torpedo. Success of the attack depends 
primarily on the availability of maintaining azimuth-plane 
steering to within a short horizontal range of the target, 
and the subsequent operation of the influence exploder. 

Case A is of particular interest.to this analysis, as 
both the search and attack-depth are within the subsurface 


layer. 
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Ift~. OCEANOGRAPHIC BACKGROUND MATERIAL 


The oceanographic background for predicting typical 
enawextreme conditions of 

-wind speed 

-wave heights 

-bubble densities and distributions 

-ambient noise versus self noise 

-sound speed profiles 
are outlined in detail in Appendix A. Even though most 
data have beneral validity for Norwegian coastal waters, 
the region above 68°N are of particular interest. fThus, 
a typical area combining open as well as confined waters 
can be represented by "Andfjord" at 70°N, where the oceano- 
graphic conditions can be related to the weather station 
"Andenes," see Fig. 2. 

Figures 3 and 4 [Ref. 1] show average windspeed and 
the occurrence of significant wave heights as a function 
of time of year at weather station "Andenes," respectively. 
The bulk of data is centered around a windspeed of 
Beaufort:4-5 (11-21 kts) and SS:3-4 (significant wave heights: 
1-2 m). Table I gives the relationship between SS, wind- 


speed and expected significant wave heights. 
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TABLE I 


RELATIONSHIP BETWEEN SEA STATE (SS), WIND SPEED, 
AND EXPECTED SIGNIFICANT WAVE HEIGHTS (He) 


Class Significant Wave Height Beaufort 
Number oie Scale 
SS He 
0 0) 0 
ih 0 - 0.1 il 
2 OnLy ="055 2 
B 055 -.1.25 3-4 
4 lie25 = 245 5 
> 2.5 - 4.0 6-7 
6 ia — wo. 0 8 
7 6.0 - 9.0 9-10 
8 9.0 - 14.0 ae] 
3S >14 eZ 
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This, together with the low probability of having an 
amnibious operation occurring in high sea states (SS>5) 
justifies the assumption of moderate sea state with wind- 
speed in the region of 12 kts and wave heights of 2 m. 

The bubble data distribution taken from Ref. 2 was 
obtained in the area "Troms6" - "Bjdérnéya" during the 
period June-November 1978. These data correlate very well 
with a larger body of data obtained by H. Medwin [Ref. 3]. 

Figures 5, 6 and 7 [Ref. 2] show the density of resonant 
bubbles as a function of depth with windspeed as parameter 
for the 12,38 and 120 kHz. As seen from these figures, 
the number of resonant bubbles are an increasing function 
of frequency and windspeed, and a decreasing function of 
depth. Below a depth of approximately 15 m, the number of 
bubbles 1s negligible for the windspeed of interest. 

The effect of SS (windforce) on the ambient noise level 
1S given in Fig. 8 [Ref. 4]. Shallow coastal Norwegian 
waters are typically 5-10 dB noisier than the corresponding 
deep water. However, great variability caused by local 
Ship traffic, fishing fleet activity, marine life and 
local wind conditions makes ambient noise prediction diffi- 
cult in these areas. This means that for accurate modeling, 
ambient noise prediction have to be done at each location 
as its level is both site and time dependent. However, 


Fig. 8 shows that for frequencies higher than 50 kHz, the 


ng 





effect of windforce on the ambient noise level decreases 
to a lower bound determined by the thermal agitation. 
Based on the above discussion and experience related to 
noise levels for torpedo systems, the self noise will be 
assumed to be dominant through this analysis. 

Figure 9, obtained from Ref. 5, shows that the sound 
speed profiles usually encountered in the area of interest 
Pesim@ts in extremely difficult sonar conditions. This is 
illustrated in Figs. 10 and 11, which show worst-case ray 
path derived from Fig. 9. In addition, the presence of 
bubbles in the subsurface layer causes the sound speed to 
be a function of frequency. The above two factors may 
frequently be the ones limiting the detection range of the 
torpedo. These effects can be minimized by selecting an 
appropriate search depth for a particular sound speed 
profile. In addition, for Case B the corresponding curved 
homing trajectories in the pitch plane give an error in 
apparent range to target. This effect will normally be 
taken into account by devising appropriate attack logic 


which is outside the scope of the present analysis. 
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IV. THE PASSIVE SONAR EQUATION 


A measure of efficiency for a passive homing torpedo 


is the detection range obtained by solving the passive 


sonar equation for broadband noise: 


where 


SL+10logB-TL (geom) -aR+DI-NL-DT=0 ) 


SL 


DI 


NL 


DT 


spectral level of the broadband noise 

radiated by the target (in dB re luPa/lHz at 
le 

detection range (in m). 

attenuation coefficient at the center frequency 
(Tn asm) 

receiving sensitivity (directivity index) 

(1n GB re luPa). 

noise level at the receiver in the bandwidth B 
(in dB re luPa). 

detection threshold; the signal to noise ratio 
at the transducer output required for a de- 
Beccien probability Of Pp and associated false 


alarm probability Pra (Gina Bye 


Zu 








V. REFERENCE MODEL 


oo. INTRODUCTION 

In order to produce the reference for the analysis 
the sonar equation is solved assuming ideal free-field 
conditions, a simple noise source model, and a generalized 


square-law detector. 


B. IDEALIZER TRANSMISSION LOSS MODEL 

Because the presence of refraction, scattering, and of 
ocean boundaries, free-field conditions associated with 
homogeneous (isovelocity) and unbounded medium seldom 
exist in the sea. However, as a basis for comparison, 
the ubigutuous spherical spreading law plus an added loss 
term due to "normal absorption" can be used as a reference 
model for measuring the effects of the previous mentioned 
scattering and absorption mechanisms. Thus, the reference 
transmission loss model can be expressed as: 

TL=20logR + aR C2) 

where the absorption coefficient, expressed in dB/m, can 


be obtained from Fig. 14 taken from Ref. 8. 
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oe MOlSE SOURCE MODEL 
1. General Characteristics of Noise Sources 

Sound 1s generated in a fluid medium by any process 
that causes an unsteady pressure field. Physically processes 
that can cause an unsteady pressure field include: 

-Pulsation of a boundary surface of the medium 

-The action of a nonsteady source on the fluid 

Zrouremt MOELON in the fluid 

-Oscillatory temperatures 
It can be shown, e.g., Ref. 10, that each source mechanism 
mathematically corresponds to a dominant order of multipole. 
If all sources are of such a nature that their time variation 
can be described by a Fourier Integral, it can be shown 


[Ref. 10] that the Helmholtz Equation is 


2 
z W or oe 
V P,, (x) ~ cr Pp, (x) = 47f£ (x) (3) 
30 (x) A = ees 
Se 
es) 
Term 1 Term 2 Term 3 


where the right hand side describes distributed source terms. 
Mae terms on the right hand side of Eq. (3) have the 
following interpretations: 
Tea |: mass injection 
Term 2: external force 


Term 3: turbulent shear stress 
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In the long distance and long wavelength approximation, 
it can be shown that the mass injection term gives rise to 
a simple source; a zero order pole called a monopole. The 
monopole radiates omnidirectional with no angular dependence. 
At large distances the pressure field from the monopole 
radition is that of a point source. Examples of this are: 
-Pulsating bubbles 
-Cavitation 
The external force, in the long distance and long wave 
length approximation, is associated with a dominant dipole 
which has a coSine directional pattern. Examples of this 
type of radiation is that caused by the vibratory motion 
of an unbaffled rigid body. 
Radiation from turbulent shear stresses is charac- 
terized by a lowest order term of quadrupole nature. 
The efficiency of the source terms decreases with 
increasing dependence on the spatial derivatives. This 
can be understood when recognizing that wave functions of 


the general form f(x-ct) have a time derivative 


|scf(x-ct) | = cf, (x-ct) (4) 


which is magnitude c (sound speed) greater than the spatial 


Gerpivative 


0 = _ 
laf (x-ct) | = F(x CE} (5) 
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Other factors being equal, the radiation from an external 
force is small compared to that from mass injection, and 
that from turbulent shear stress is the smallest; therefore, 
moncepole radiation is the dominant term. 

Propeller cavitation, when it occurs, is usually 
the dominant noise source for any marine vessel. Submarine 
and torpedoes often operate at a depth great enough to avoid 
cavitation. Surface ships, on the other hand, generally 
have well developed propeller cavitation with the result 
that their radiated spectrum from 5 Hz to 100 kHz is 


controlled by this source. 


The basic phenomena of cavitation combined with 
propeller hydrodynamics give the fundamental characteristics 
of propeller cavitation noise. An excellent qualitatively 
discussion of this can be found in Ref. 9:Chs. 7 and 8, 
from which the following is extracted: 

Propeller blades are rotating twisted wings 
that produce hydrodynamic forces. Depending on operating 
conditions, they experience cavitation on a number of 
different places. Of these there are three prominent types: 

-Tip vortex cavitation 

-Hub vortex cavitation 

-Surface blade cavitation 

In addition to the two types of vortex cavitation, 


there normally are two types of blade surface cavitation: 


Zo 





mecdex; driving face 

Zoeome: SUCTION surface 
Of all kinds of propeller cavitation, surface blade cavi- 
tation on the suction surface is normally the most noisy, 
While hub vortex cavitation is the least noisy. 

2. The Noise Source Model 

Due to lack of recorded and available noise data 
from the target in question, the noise source has been 
generalized on the basis of the following discussion and 
assumptions: 

The noise source will be build up around a surface 
blade cavitating propeller operating in a good to poor 
wake; surface cavitation will be assumed to be dominant. 

D. Ross [Ref. 9] has developed an approximate 
theory for cavitation noise, where dimensional analysis 
1s combined with the basic results of cavitation theory 
that the acoustic pressure 1S proportional to the product 
of the collapse pressure of the cavities and the volume 
of cavitation produced per unit time. From this synthesis 
it is found that the total acoustic intensity varies as: 


eos e (Ui e. oU) U 
I= K,, 2S (i -)7) (6) 
ee ee era 


where 


distance of the hydrophone from the source 


ry 
I 


Key = the cavitation inception parameter 
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UL = blade tip speed 


UL; = blade tip speed for inception of cavitation 


This expression shows that propeller cavitation noise power 
1s proportional to the total number of blades, b, the blade 
chord, s, and to the propeller diameter, D, and is a function 
of the tip speed with the dependence on the tip speed being 
the strongest. The different blade surface sections where 
cavitation exist are uncorrelated and the radiated noise 

1s treated as a single monopole radiation so that ata 
distance r>>a (where a is the characteristic dimension on 

the source region) the radiation is similar to that of a 
point source with no angular dependence. 

Submarines and torpedoes with centerline propellers 
have a relatively symmetric inflow condition. Surface ship 
propellers, in contrast, operate under nonuniform inflow 
conditions. Circumferential wake variation causes the 
radiated sound to be amplitude modulated at the blade rate 
frequency. Furthermore, slight physical difference 
between the blades produces modulation at the shaft rate 
frequency. These amplitude variations gives a very 
distinct characteristic to the radiated noise that can be 
used for classification purposes to reduce the probability 
of false alarm. 

The most complete source of data on surface ship 


radiated noise are measurements made during WWII, reported 


2G ) 





iiecaecompendium issued by the U.S. Office of Scientific 
Research and Development (OSRD) in 1945 and declassified 
in 1960. When these data are examined the radiated noise 
is found to depend on tip speed and the number of propeller 
blades with little dependence on the other variables. For 
surface ships near cruise speed, the source level for 
frequency over 100 Hz can be written as: 

Sia '+Z0=-20llogt; £>100 Hz C70) 
where 

f = frequency in Hz 

SL'= overall source level in dB re 1l1uPa. 


The overall level can be expressed as: 


SL'=175+60logU /25 +10logb/4 (8) 
where 

U = mnD (9) 

n = rotational speed (rpm) 


D = diameter of the propeller (m). 
The above expressions are used as the basis for the noise 


model with the following input data: 


n = 300 rpm for maximum cruise speed of 15 kts. 
n = 180 rpm for a cruise speed of 10 kts. 

ea—ee m- 

ba = >. 


The resulting noise spectrum, in dB re luPa at lm, 


as a function of speed in kts. are tabulated in Table fui. 
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and plotted in Fig. 15. A one sigma region (5 @B of un- 
certainty) is incorporated in Fig. 15. 

The above noise spectrum estimation agrees very 
well in the high frequency limit, with more recent studies 
by A. Lovik [Refs. 1l and 12]. Here it is found that the 
cavitation spectra, both theoretically and experimentally, 
can be divided into four frequency regions, as illustrated 
ma Fag. 16. 

Region I is dominated by noise at the blade fre- 
quency and its harmonics. The emitted sound is caused by 
the volume variation of the main cavity. 

Region II starts at the bubble frequency, which is 
the reciprocal of the lifetime of the main cavity. The 
mean power level is found to decrease with increasing fre- 
quency as cee. 

Region III is an intermediate region. 

Region IV associated with the shock waves starts 
at the mean collapse frequency Ear given by the mean 
collapse time. The power level is found to decrease as 
€*, as in Hey. 5. 

The number of gas bubbles in the water have a 
pronounced effect on the high frequency cavitation noise 
Mmeom tae propeller. This is illustrated in Fig. li7, 
obtained from Ref. 12, where the power is found to decrease 


as much as 40 dB with increasing gas content. 
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scaling laws are developed [Ref. 12] for each region 
based on a series of models and full scale measurements. 
These laws depend on the dynamic pressure induced by the 
propeller, the model ratio, and the gas content of the water. 
The full scale measurements were performed in cooperation 
with the Royal Norwegian Navy and the Marine Institute of 
Norway. The model experiments were performed in the largest 
cavitation tunnel at the Ship Research Institute of Norway. 

In summary, the scaling of cavitation noise was 
demonstrated to be a useful tool in predicting a full scale 
cavitation noise as shown in Fig. 18 [Ref. 12] which compares 
measured noise spectra for the model and full scale 
measurements. 

For the high frequency region, the source levels 
are of the same magnitude as predicted by the WWII empirical 


formula. 


D. PASSIVE MODE RECEIVER CHARACTERISTICS 
1. Assumptions 

For receiver characteristics assume a square law 
detector with a center frequency f£-60 or 30 kHz and a band-~ 
width B. The detection scheme is shown in Fig. 16. The 
principal assumptions employed in the derivation are as 
follows: 

~Gaussian signals in Gaussian noise 

-~Frequency independent signal and noise spectra 


ou 





-Integration time T is sufficiently long to 
permit application of the central limit theorem. 
Ze ber ivations 
MeowcGeuSector input r(t) is assumed to be a zero-= 
mean Gaussian process composed of noise alone or signal 


plus noise expressed by the two well known hypotheses 


Beit) = n(t) (EO) 
H,yr(t) = s(t)+n(t) 
where: 
n(t): noise signal 
Se): signa. 


The two signals s(t) and n(t) are assumed to be independent. 
Assume that the spectral shape of s(t) and n(t) 
are the same, such that Hy and Hy Onis vatirecr in ene, cotal 


power level. Then the detector-smoother have the form: 





Schematic of detector-smoother 


and 
eae) =) r- (t) (11) 
Furthermore, let the noise variance be normalized to unity 
(for convenience) and the signal variance by denoted by o 
Var[n(t) J=1 (a) 


Merie(t) |=o- 


a 





Because of the assumed similarity in the spectral shapes, 
the autocorrelation functions are 

R(t)=F[N(£) J=9(t) ; N(£) is the noise power (13) 

Spectral density. 

Paw) okt) 

ao) H.e(t) 

Hy (l+0*) 9(t) 

Furthermore, assume that the integration time T is long 
enough so the central limit theorem holds, implying that 
K also is a Gaussian random variable. 

This yields that the probability density function 
of the output variable and hence the detection and false 
alarm probabilities are completely determined once the 


mean and the variance of K are derived. 


If a process V(t) 1S wide-sense stationary, then 
a 
E[v(t) ]= of (atte poe kh cae ale (14) 


Thus, assuming that r(t) is a wide-sense stationary process. 
E[K]=E[x(t) ]=E[r(t) ]=1l+0° ins) 


and similarly 


a 
var[v]=E(v?]-<Elv]}?= 4, // Blv(t)v(s) ]atas-v? (16) 


cl TE 
iE 2 al a 
Coe =F [/ [R.(t-s) -v Jdtds=—; | Cov. (t s)dtds 
00 00 
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where Cov, (t-s) 1S a covariance function. 
Then letting 


T = t-s 


C tts, and substituting into Eq. (7) yields 


5 
2T 
var [v]=3, //cov, (1) S58 ; 
oO 7 T TF 
a 
var[VEz / (i-th icov, (ry at (17) 
ee 
Consequently 
al 
var [K]=2 J a-tphicov, (nat (18) 


=i 
Then, evaluating the covariance function from the auto- 
Gerrelation function 
R,(T)=E[x(t)x(t-t)] = E[x* (t)r?(t-t1) J (19) 
Since r(t) is Gaussian, the above fourth moment can be 


expressed as product and sums of second moments: 





es (3) er7 (ttt) + 2r(t)r(t-t) + r(t)xr (tT) 
oe 2 
ea) ~ 2R_ (tT) 


Bet) (1+ °)°+2(1+0°)° Oo tr) (20) 
Thus, the covariance function is 
a = iz 
Go (1)=R, (7) fis seal} 


=(ito-) “te l+o-)o-(t)—(l+o-) 
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Sac) =ea@iera*) *p* (7) (21 


Inserting Eq. (12) into Eq. (9) yields 


dl 
var(kjJ=e / ¢(1-lehy 7a N( dLtpor))| ear (Ge el eke 
onl 


Ee 
Zy 2 
var [k]=20t0 ) f (i-tely p*(t)dt (22) 
-T 


If T 1s large compared to the correlation time 
WB > > 1 


ea@ewewe Can substitute the limit for Eq. (22) by 


oO 


zy) 2 
Var[K]= fio /[ pO acral (23) 
OHO). [ snes 
aT 


It we further make the assumption that the signal and noise 
have ideal flat bandpass spectra: 
precy = | 1/2B, £=-B/2<£<£+b/2 (24) 
otherwise, 


Pmieetteamg Eq. (24) in Eq. (23) yields 


Cosy 


=a Zi) 


Var [K]= 


The probability density function for the output of the 


detector have the following form 
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=Threshold 


Ho choosen aa eS H, choosen 


Here w is the outcome of all possible signals. The false 
alarm probability is obtained by integrating the conditional 
probability p(w|H.) over the outcome space for which to 


choose Ht, 


ie J 
Pra = i p(w] H,) dw = Q( 
5 





) (26) 


Similarity, the miss probability=l-detection probability 
is obtained by integrating the conditional probability 
p(w|H,) over the outcome space for which to choose Ho: 


J 
u 
Py = 1-P, = | p(w|H,) dw = Q( ee ) (ZI) 





Further defining the input and output signal-to-noise 
ratios as 


S/N(input) = 2 = o? (28) 


S/N(output) = ———— (29) 
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For the square law detector where 


uy = E(n*(t)] =1, u, = Blr?(t)] =E[K] =1+02 (30) 
oe = Var[n*(t)] = aa OS ewer (Ge) | 

7 me (koe) 

= Var[K] = 37 


The output Signal-to-noise ratio: 


2a 4 
S/N (output) = EC aah) = = BT oO (aa) 


a (1+0°) * 
The input signal-to-noise ratio: 

S/ Monee) = o7 (32) 
The probability of false alarm: 

Bex = Ql yBr(j-1)] (33) 
The probability of detection: 


ane os 
py = 1-o[vBT (=25))] (34) 
I+o 


Mme E@uations (32), (33), and (34) are plotted for a 
variety of realistically encountered conditions. Figure 17 
gives BT versus S/N(in) @dentical to the detection 
threshold DT) for various combinations of Pp and Pra: 
Figure 18 gives Py versus S/N (in) Sor yas lous 
combinations of the threshold j and Pra: This constitutes 


the ROC-curves for the generalized square law detector. 
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E. REFERENCE DETECTION PERFORMANCE 
i “lncrodmetion 
As both the scattering and absorption are frequency 


dependent, it is necessary to have reference models for both 


60 and 30 kHz. 
2. 60 kHz Case 
The range dependent portion of the sonar Equation (l): 
-20logr=-GR 
MeepLOtted in Fig. 19 for a=0.021 dB/m taken from Fig. 14. 
Figures 15, 17 and 19 are then used to estimate the 
detection range: 


(a) Select Pp: B, and the integration time T. 


Poa! 
The detection threshold DT=S/N(input) is then 
feud irom Fig. 17. 

(b) Select the speed of the target and find the SL 
from the noise source model (Fig. 15). Then, 
reasonable values for the receiver sensitivity 
DI and the self noise level NL yields the left 
hand side of the sonar equation (1) except for 
range dependent term. 

(c) Use Fig. 19 to solve the passive equation 
oO rs 


A realistic example may illustrate the above procedure. 
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(a) Entering Fig. 17 with 


Ph = OFFS 

PRA ~ 107° 

B = 4500 Hz 
T = 100 msec 


yields 
pi = =625 dB. 
(b) Selecting a target speed of 12 kts. gives 
(from Fig. 15) SL=100 dB. Selecting a 
typical transducer sensitivity DI = -180 dB. 
Assuming the NL to be dominated by self noise 
of typical value NL=-124 dB. This yields 
SL+10logB+DI-NL-DT=87 dB. 
we) “P2oqure 19 with 
-20logR-0.02R=-87 
yields 
R=1200 m. for a=0.02 dB/m and £-50 kHz. 
The influence of different design parameters like self 
noise and detection threshold on the passive detection 
performance is now easily investigated by the above 
procedure. 
Although outside the main scope of this analysis, 
the above statement can be confirmed with an example. As 
seen from Fig. 18, a probability of false alarm p,, = om 


implies a threshold setting = 0.9 dB above the noise level. 


a 





Due to the variability of the noise level together with 
@ae practical difficulty in accurately setting the threshold, 
a more realistic goal for the threshold would typically be 
3 dB. Going into Fig. 17 shows that the corresponding value 
Bor DT for Ph = 0.5 is DT = 0 dB, with a corresponding low 
value for the Pra: 

Letting DT=0 dB and keeping the previous assumed 
values of SL, DI, and NL yields: 

SL+101l0gB+DI-NL-DT=80.5 dB. 

The corresponding detection range is: 

R = 1000 m, for a = 0.02 dB/m and £ = 60 KHz. 
Thus, this change in threshold setting caused a decrease 
in detection range from 1200 m to 1000 min return of a 
Significant decrease in the false alarm probability. 

S. 20.kHz Case 

In order to estimate the reference detection range 
for an operating frequency of 30 kHz, we utilize the sonar 
equation (1). 

Assume that the receiver has the same generalized 
passive detector characteristics as in the 60 kHz case: 

pia —180. dB 

NL = -124 dB. 

DT = -6.5 dB, based on Pp = 0.5 and Pea = ere. 
However, the empirical equation (7) for the source level: 


SL=SL'+20-20logf 
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Shows that the source level falls off as Ee If the 


dynamical and dimensional parameters of the propeller are 
the same, SL will increase by +6 dB when the frequency is 
reduced from 60 to 30 kHz. 

A source level of 

SL=100+5=106 dB 
gives a range dependent solution of the sonar equation 

-20logR -aR=-87-6 = -93 dB. 

AeoLOc-Of 

-201logR-aR 
1S given in Fig. 20 for an absorption coefficient a=0.01 dB/m 
Peeen Crom Fig. 14. 

Figure 20 then gives a detection range of 

R = 2400 m. 
Thus, as seen from these ideal reference calculations, 
halving the frequency gives a higher source level anda 
lower absorption loss, resulting in a doubling of the 


detection range. 
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VI. THE EFFECT OF SURFACE SCATTERING 


A. OCEANOGRAPHIC DESCRIPTION OF THE SEA SURFACE 

The roughness of the sea surface is the essence of 
mame Scattering mechanism. Thus, to adequately describe 
the scattering of sound from a randomly rough sea surface, 
aC 1S necessary to formulate a suitable description of the 
sea surface from an acoustical propagation point of view. 

Generally the shape of the rough sea surface is most 
appropriately described in terms of time and spatial 
dependent random variables. However, observation of the 
ocean under the same environmental (meteorological) condi- 
tions indicates that the roughness is the same over large 
areas and for periods of at least several hours. The 
random processes responsible for the structure of the sea 
surface, therefore, can be considered stationary at least 
over periods of hours. With this assumption, the sea 
surface can be described in terms of the statistical 
description of the surface displacement function, and the 
distribution of signals reflected from the sea surface 
can then be related to this probability distribution 
Zune tion. 

The most significant statistical parameters describ- 
ing the scattering mechanism from the randomly rough 


surface are: 
42 





-the mean square slope 
-the mean square surface height 
-the correlation length. 
Optical measurements made at sea by C. Cox and W. Munk 
[Ref. 13] showed that the sea surface with an arbitrary 
wide continuous spectrum of waves is characterized by a 
Gaussian distributed surface slope. The mean square 
slope, determined from these optical measurements 1S 
ame = = (345-120) x100 (35) 
where 
W=wind speed in m/s measured 41 ft (12.5 m) above 
the sea surface. 
The Gaussian distribution of surface slopes implies that 
the surface displacement function can be described by a 
Gaussian probability density function with zero mean 
<tT> = 0 
and varlance 
aa - 
men Gaussian correlation function 
amen) EX tr) > 
The mean square height 0”, is obtained by integrating the 
frequency spectrum of the fully developed sea. The fre- 


quency spectrum G is given by the commonly accepted semi- 


empirical expression of Piersom-Moskowitz [Ref. 14] as 
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Sia Q 
G(Q) = ay exp[-8 1" (36) 


where 


Q= frequency (in s7) 


) 
G— €. 1 

B= .74 

a : -l 

a ewe(ans ~)s 

W=wind speed in m/s at 19.5m above the sea surface 
See Gravitational acceleratiion (in m/s-). 


This gives 


4 
2 - / Gide] — aS fou.) 


0 48q° 


For a Gaussian autocorrelation function expressed as 


52577 one 
<t(t)t(t+t)>= e 1/7 F (38) 


Jr 


W(t) = 


Q 
i) 


where 
te—. correlation length. 
The following relationship for the mean square height holds 


for sea of small roughness 


ae eG | 
te = , see later Eq. (60) (39) 


me 





and the correlation length is thus: 


T= 72 = (40) 
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Pee OCATTERING THEORY 


mie tbeal boundaries are rough for radiation with short 
enough wavelength, and the apparent roughness depends on 
the "viewing" conditions. The wave reflected by a plane 
surface has the same properties as the incident wave since 
the radiation is scattered coherently and there is a definite 
relation between the incident and scattered waves. 

A randomly rough surface, however, such as the wind 
generated ocean surface, scatters radiation in all directions, 
i.e., an illuminated area is visible from any direction. 

Heuristically there are two distinct approaches to this 
phenomena. 

1. If the boundary is rough most of the radiation 1s 
scattered and there is little transmission in the 
specular direction. Thus, the attenuation caused 
by the irregularities can be included in the trans- 
mission equation. 

2. If the surface is truly smooth, it can be assumed 
that the effect of the boundary is to supplement the 
original pressure field by an out-of-phase image 
contribution. For a randomly rough surface the 
reflected sound neither completely cancels the 
direct sound nor adds to give +6 dB pressure peaks 
of the interference pattern. For a rough surface, 

this supplement is a small fraction of the direct 


path. 
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The second approach will be used with the Simplified 
assumption that the sea below the surface has an iSOLCrOpLC 
statistical description; i.e., the mean acoustic velocity 
and the mean density are assumed to be constant and have 
negligibly small mean square fluctuations. 

The estimation of the scattering is based on an approxi- 
mation method employing both ray and wave theory. Ray 
methods are used to follow the acoustic signal from the 
Meise SOurce to the vicinity of the sea surface. fThen, 
wave theory is used to calculate the scattering process. 
Finally, ray theory is used to follow the scattered signal 
to the receiver. 

The geometry is given in Fig. 24a. The origin of the 
coordinate system is at the center of the illuminated area. 
The x-y plane coincides with the mean of the rough surface 
as averaged over the illuminated area. 

The source and receiver are at distances Ry and Rot 
respectively, from the origin. Ry 1s the xz-plane and 
makes the angle oy with the z-axis. Ro makes the angle 0. 
with the z~axis and the projection of Ro on the xy=plane 
has an angle 8 relative to the x-axis. 

For high frequencies Ry and Ry are much larger than 
the acoustic wavelength. Then both the incident wave 


and scattered waves can be treated as nearly plane waves. 
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The formulation of the scattering problem will be 
based on the Helmholtz integral which requires known 
values of the normal derivatives of the incident and 
reflected waves on the boundary. These are estimated 
by using the neuristic Kirchhoff's approximation, which 
assumes that the wave is locally reflected by a plane 
surface; 1.e., an approximation restricted to a surface 
not too rough and not shadowed. 

Further, the receiver derictivity, as indicated by 
Fig. 24b, will be used to limit the surface area that is 
illuminated. 

Since this procedure is based on a detailed develop- 
Ment by I. Tolstoy and C. S. Clay [Ref. 15!, only the 
Main points will be outlined here to bring out the assump- 
tions made and the inherent limitations of this approach. 

The development starts by considering the inhomogeneous 


wave equation of the general form: 


—. 1 3°p (x,t) . 
ee) SOS EX, 8) (41) 
Sgr. toes 
where 
f (x,t) —Slcmamihownims@unrce idaistri bution. 


The development is based on the following initial assumptions: 
-The medium is homogeneous. 


-The medium is bounded by some surface S, onto which 


an incident wave impinges. 
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-The boundary is characterized by the Specific 
acoustic admittance and the shape of the boundary. 
-The incident wave is harmonic. 
The assumed harmonic source implies that £ (x,t) can be 
decomposed into a Fourier integral. Furthermore, assuming 
that the solution of Eq. (41) can be decomposed in time, 


we arrive at the Helmholtz equation 
V*p (x) + => p (x) = -47£, (x) (42) 


BEponould be noted here that L. Fortuin, in [Ref. 16] 
Showed that the Helmholtz equation is not exactly correct 
for a medium with a time dependent boundary. The equation 
can, however, be used with a good approximation when the 
time derivative of the surface elevation is much smaller 
than the speed of the waves through the medium. For under- 
water sound waves scattered by the rough sea surface, this 
means that the wind speed has to be much less than the 
sound speed; a requirement easily fulfilled for our inves- 
tigation. 

Green's method allows the solution of this linear in- 
homogeneous wave equation to be expressed in the 


heuristic Helmholtz integral form: 


i op, (x) 
pc [ £,, (2°) 6, (ax!) at tgp (6 2x0) 2 (43) 
Vv ee n 
8G ( xx") 
- ae 
n 
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The first integral on the RHS of Eq. (43) contains 
the sound sources and the bulk (volume) scattering. The 
second integral represents the surface scattering and is 
taken over all finite surfaces. 

Now, disregarding the direct path, the signal at the 


receiver is given by the surface integral alone: 


dp 
\S) Berk: l Ww 1, 0G 
Pi es) -7— {G(x2x") int = P(x ar da’ (44) 


Z on on 


where the subscript s denotes the scattered field. In 
order to solve Eq. (44) the following must be done: 
-~Give an approximate expression for the incident 
wave. 
-~Find an appropriate Greens function. 
-~Make an approximation for p(x") and dp, /dn at 
the surface. 
As we already have assumed a simple harmonic source, 
the incident wave can be expressed as: 


(i) 
Py 


(x1) = (ABS) * B etKR = BD QtkR (45) 
where 

I: power output. 

Perl lamination function. 

Assuming kR>>l, i.e., that the distance of the source 


is large compared to the wavelength, the wave in the 


bounded ensonified area can be considered as a plane wave 
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characterized by its propagation vector 


—> 
xX. 





“a+ 


thus obtaining the expression for the incident wave: 


en BD ikR; iki -x! 


P aC e (46) 


u) 
Further assuming that the receiver is sufficiently far 
removed from the scattering area with the rest of the sea 
surface not contributing, then the scattering area acts as 
a small induced source in a free space and we can approxi- 
mate the propagation of the scattered waves from the en- 
sonified region in terms of the free field Greens function: 


1KR > — 
Bes l 


2 





G(xX2x1)* (47) 


where 7 
x, 
k = k-= 
Ss ee 
We further assume that each surface element da’ acts as a 
small reflector, and that the response of da’ to the inci- 


dent wave is that of a "local reaction," i.e., independent 
of any other part of the ensonified area A. 

Then, the Kirchhoff's approximation where it is assumed 
that p and 3p/3n vanish everywhere on the surface except at 


the ensonified area and that the values of p and dp/dn are 


proportional to the incident wave, allow the scattered 
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"target" strength for the ensonified area to be approxi- 


mated as: 
p(x )| = ce) 4 
Wy i R Py = (48) 
a 
dn? A 3n! 
where 


p is the locally reflected wave. 


seeleetion coefficient. 


1 


_ p'ctcos @-pce cos 8! ¢ a ee 
= ep e Cos Y=pCc | ' = 


sind 
p'c!cos 9@+pc cos 8? sing! 








Finally, assuming a Gaussian illumination function: 


x? vie 


2 
eso !!!C (49) 
where 
X and Y are the effective dimensions of the illuminated 


abcd, 


The scattering integral based on the Kirchhoff's approxi- 


mation can be expressed as: 


1k (Ri+R2) 
(s) _ _ikBe ri 2i (ax? +By? ) 
er (x2) = TRS — RE (8 1 8283) Jr 
atid (XY) aytayt (50) 
where 
Qa = < (Sin 6, - Sin 92 cos §3) 
g =-5 Sin 8» sin 83 
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a -S(cos oy COS G5) 


c(x'y') = surface displacement function. 


The above scattering integral is then applied to a randomly 
rough surface where the surface displacement function ¢ is 
a random variable assumed to be represented by a Gaussian 
PDF expressed as 


a: 2 
i at 2g #29 S)2L)) 


oV2T 
Math zero mean and variance o?*. 

Also, assume that the surface is slowly varying so that 
the signal reflects from an essentially stationary surface 
and that the succession of received scattered signals Ph 
are assumed to form a satistically independent set of 
Sample functions from which sequence N the first and second 
moment of the field can be investigated. Doing so Clay and 
Tolstoy [Ref. 15] found that the mean reflected signal can 


be expressed as 


-2 ie f 
<p, 5) (x2) > = Eee ¥ . (52) 
where 
_ _(s) 
Po ~ Py (x2) | a0 


is the signal reflected by a mirror-like surface, 


O©bmer factors being the same. 
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It 1s seen from the above that 

1. For o-0, all displacements cz have zero proba- 
bility and the mean signal tends to Po: Further- 
more, all elements contribute to the scattering 
coherently. 

2. For o>>k, all displacements z are equally probable. 
There are large phase shifts between contributions 
from different surface elements and they tend to cancel 
each other and the scattering radiation is incoherent. 


The second moment is defined as: 


: o pee 
<s*(t) S = va See ee <PP*> =< Dp (53) 
6 le 
where 
ee 
a 2 a 
2° ye) lt P(t) N<p(t) >] (54) 


In this expression the operation of squaring the signal has 
to be considered. However, each surface element has a 
different ¢ for a random surface. Thus, the probability 

of finding element dx' dy' with 71, and element f2 with 

dx" dy" is expressed in the bivariate distribution function 


assumed to be GausSian and of the form: 


W (102) =—— =e [-—= (6174017-201020)] (55) 
2707 (1-7) 2(1l-*)o? 
where 
v(E yn) = = <t,(xty!t) Cy (x"y"t)> (56) 
Oo 
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1s the cross-correlation function characteri Zing 


the surface shape. 


Changing to polar coordinate leads to he Oso. OL Retr. 25 


where: 


oO 


sya. 2 Re 
ee*ea f po (2er)[e4Y % (AW) 49") vay (57) 


Here J, 1s the Bessel function of zero order and x is the 


transformation parameter given as 


(or 
cos 86 





As seen from Eq. (57) 

1. For rough surface, y*o* is large and the second 
term in the bracket, the coherent part, is 
negligible. 

2. For smooth surface, y*o* is zero and the whole 
bracket 1S zero. 

Since Eq. (57) cannot be integrated directly, Clay and 
Tolstoy [Ref. 15] consider it for small and large yo 
Separately. 

Taws, concentrating on the high frequency limit, Clay 

and Tolstoy assumed y*o*>>1l and showed that: 

-The coherent component is negligible. 


-<s*> =<pp*>, as the means tends to zero, 
oc 


és? >S<ppt>a fo, (2x 
0 


222 SS 
None Ome eee ln ( 58) 
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Equation (58) consists of the product of an oscillatory 
function and an exponential function. Because of the Bessel 
function, the main contribution to the integral is near r=0. 
Near r=0 the phase changes slowly and the expression can be 
evaluated by the method of stationary phase. Thus, the 


expression of ~ about r=0 is given as: 


Zz 
b x 1+y" (0) 5 (59) 


Furthermore, Clay and Tolstoy show that ~"(0) can be related 


to the characteristics of the surface as: 


}y"(0)| = = <t'*> (60) 
G 


Finally, Clay and Tolstoy show that the scattering Signal 


can be expressed as: 





AG se od 
<s*> = <pp*> = <P, BR*> Sa; (61) 
ee 1 
eee 
where 
A = ensonified area 
<p,p#>: the expected average value of p,*, where 


p, is the incoming pressure to the 


luminated area. 





Sue = scattering function 
x? 
2 2 = 
_ bc) 


Diteos 0, + COS On mG me 
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1+cos8,cos0.-sin@;sin@2cos6; 


ey e282 83) = cosé,+cosé@ 
1 we 


1 


ee 
p . 0 © cos® - pc cosé'. coe. ec 
; 2 Were 


p'c'cos8+pc cose! sing sing? 
eee, 
Ve ai ni (cos @,+ cos @>) 
a 2 
cos @ 





Qa = = (Sth 9,- Sin 8. cos @3) 


2 
we a2 
i r 
<g'*> = (3 + 5.12W) x ome W = wind speed in m/s. 


Bemaeion (61) is valid for: 
ayeG >>I 


ayt<ct?sy>> de 


In summary, the reflection of high frequency signals yields 
Scattered radiation which is incoherent. Furthermore, as 
pointed out by Clay and Tolstoy in [Ref. 15] although the 
radiation is primarily scattered in eee eee direece Tony 
parts are scattered in all directions. As seen from 

ma. (62) the scattering function Sue 1S primarily dependent 
on the mean square slope of the surface <c’*> and neither 
the mean square wave height o* nor the correlation distance 


are important. Furthermore, it is noticeable that in the 
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high frequency limit the scattering function is independent 


of the frequency since 


e- a +1 when r7>+0 
Cee) 

C. GEOMETRICAL SHADOWING 

The phenomena of shadowing of some surface areas by 
others has to be considered either when the surface irreg- 
ularities are large with respect to the incident wavelength 
or when the grazing angle is small. 

The few papers devoted to this subject are aimed mainly 
towards calculation of a "Shadowing function" based on 
the statistics of the surface. 

An explicit method, geometrical shadowing, has been 
introduced by P. Bechmann [Ref. 17] where the shadowing 
function $(8) is the probability that the point ¢ (Fig. 25) 


jms illuminated. 
S(§) = exp t= f a(x ax) (62) 
0 


where: 
q(x) is the probability that ¢ is shaded by ¢ 27) 
the interval (x,x+tdx) given that it is not shaded 
Beran) . 

This calculation of $(8) only considers the elevation of 

the surface observation point. However, the slope also 


plays a role in that if its value exceeds cotd the point 
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will certainly be shaded. Thus, R. Wagner, in Ref. 18, 


incorporated both ¢ and ¢' using P. Bechmann's method and 
found that the conditional probability that a point on the 
Surface 1s illuminated, given that it has height ¢ and 


i 
Slope ¢', can be expressed as: 


en e2751°) = exp[- fa (x) ax] u(cot@ - z') (63) 
0 


where 
u: is the unit step function. 
eer: 1S the cenditiomal probability that &@ is 
shadowed in the interval (x,x+dx) given that 
it is not shadowed in (0,x). 
The function q(x) cannot be calculated exactly. Thus, 
R. Wagner made the approximation that, for all x, the 
probability that 7 crosses the ray in dx is independent of 
Bae Values of t and c* at x=0. In the above, no mention 
has been made of the direction of obServation. However, 
in this respect, R. Wagner [Ref. 18] points out that in 
the high frequency limit only those portions of the 
surface which are illuminated simultaneously by rays in 
the direction of incidence and observation can contribute 
to the observed scattered power. For this condition, 
R. Wagner introduced both directions as independent 
variables in the so-called bistatic shadowing probabilities 


where he defines: 
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alg, SMe OiG 17 Sa ' ) GO be Ene conditional propabs 11 ty 
that the surface will not cross the incoming ray 
(Ray 1) or the outgoing ray (Ray 2) anywhere, 
given that both rays pass through an arbitrary 
point or the surface with displacement t and slope ¢’. 
Qewee co; 192,01,061') to be the conditional probability 
that the surface does not cross Ray 1, given that 
1t does not cross Ray 2 and that both rays pass 
Emeough the point ¢ having slope z'. 
Thus, the conditional shadowing function can be expressed 
as: 
Biron Si,Cr') = S(8:(@2,61,61") Si(82ie1,61°) (64) 
The shadowing function is then obtained by averaging 


over all possible heights and slopes 


S(@1,82) = ff scor,ezle.c% W(t,5')de az' (65) 
0 


Memes w(z,c') is the bivariate PDF of the surface height 


and slope, assumed to be Gaussian 


rac, Gc.) = = volvo") exp{-f— “atoT (66) 
where wy, =0% and ¥'' are the values at t=0 of the correla- 
tion function and its second derivatives, respectively. 

For the region 0<O<5 where the probability of crossing 
one ray is assumed independent of that of crossing the 


other, R. Wagner found that the bistatic shadowing function 


ao 





could be expressed as: 


{l-exp[-2(B,+B.)]}xf{erf v,terf v2} 





S(9;,;,0 ee 
(61,902) T(B14B>) (67) 
where 
exp(-v,°)-/T v. erfe v, 
Bs = = > 1te=41,2 
ait 
1 
ne In | 
=<. IS 2 
207] 00" | 2x2 


as we from Eq. (60) have that o7/y"(0) | 


= ao > = re 
i = Cee 0: 12 = 1,2 
and noting that 
| 2 eat? 
Error function erf(x) = — Ve e ehe 
VT 
| 2 
Complementary error function erfc(x) = = a e Gite 
x 


The shadowing function S$(8,,92) is, in short, the fraction 
of the surface still illuminated. As seen from Eq. (60), 
the scattered field, in the high frequency case, 1S pro- 
portional to the illuminated area. Hence, the shadowing 
effect of a rough surface can be introduced by multiplying 


the ensonified area A by the shadowing function S(8,,82). 
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D. ESTIMATING THE SURFACE SCATTERING EFFECT 

The following estimations are based on calculations 
in the specular direction, which is, as pointed out 
earlier, expected to give the maximum supplementary 
scattering effect. Hence, in the specular direction where 
8, = 35 = 6 and 683 = Ome the scattering function Sue 


reduces to the following expression: 


£* (6) R? S fSeeec 
Sue SSeS D2 <7! 25) (68) 
2 (cos6,+cos6,)7<c!*> 
pe 
S = 
me peels 
as 

£(3) = cosé 

Y = kcosé@ 

a = 0 

— o — 
* ~ cose s 


As pointed out in the previous paragraph, the reflection 
of very high frequency signals by the sea surface yields 
scattered radiation that is incoherent under the assumption 
that 

eae >> 

(2y2<c2>) ood 

R 
Before we launch into the calculations, we will verify 

these criteria for the frequency range of interest: 60 and 


30 4sHZ. 
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Ue zang 





y? = k? cos7@ ao = 8.1 x 10° 
2. aw" 
= ae B = 0.74 
g 
ee = (aS -3 2 
r .12W) x 10 g = 9.81 m/s 


and assuming a windspeed of 10 m/s (SS3) ylelds 


fe = ©2350 m* > o = 0.55 m 


eee asa x 10°" 


For 60 kHz(A} =.025 m) and low grazing angles, e.g., 9 = 85° 


4y?7o2 = 5.75 x 10° >>] 


re <c'=>) = 5.2) x 10° >> + 
R?2 
and the first order | R20.5m. 
R¢ 


ie, tom the 60 kHz case, the criteria are fulfilled. 


For 30 kHz (\=.05m) and 9=85°. 
2 


ae = Ade 10°55) 
il 
(2y?<7**>) = 13>> oF 


Meaem the first order R>1.0 m. 
Thus, also for the 30 kHz case the criteria are fulfilled. 
Similarly, for the shadowing function in the specular 
direction 
O<8@,; = 62 = @ < T/2 and 83 = 0 


we obtain the following simplified expression 


G2 





$(8) = $(@,,9>) = iicesp(-4B) Jeriv (69) 


as 


ey Se Lin 


(207 |¥ "| % aac 
fe 2 = hy = Cote 
me cy =p = lexpisv")avn_v erte vj 
4/n v 


In summary, for specular scattering at the high frequency 
limit, the expected average value of p*, where p is the 
pressure field at the receiver, is then obtained from the 


following simplified expressions: 





o> )= <pp*> Eos (70) 
2 ft de 
R2 
2 
<s*> = <De Ne ae 
SuR. <c)4> 
where 
_ l-exp(-4B)erfv 
= ea 
(=e 
wang = COLO 


ae [exp (-v*) -/mv erfc v] 
me V 
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=6 > = 2° = (3 + 5.12W) x non ; W = windspeed in 


m/s. 
iat 
Pee OmcmCOSo 6] 9 GCOse: 6. = Je: 





‘ 


9'elcoss + 9 c cos9! sing 


sing! 

From Ref. 19 we use the following air/sea water inter- 
face data: 

mens aieeec’' = 343 m/s and 9 c’ = 415 Rayls. 

2. For sea water c = 1500 m/s and oc = 1.54 x he Rayls. 

The estimation of the illuminated area A for specular 
scattering where 9, = 8, = 0 and 63; = 0, is based on the 
geometry illustrated in Fig. 24b. The illuminated area 
1s given by 

A = tab (yall) 


Assuming the following data to be known 


Ae = depth of the source 
h = depth of the receiver 
a 
Ag = the half beam width of the directional receiver. 


Betn a and b of Eq. (71) can be calculated in terms of the 


detection range Ry as follows: 


h 
ole © (aes 


ee Behe 
+h 
s 2x 


Rp 


— r ° = ae Dae 
x1 = Ry - x2 = Rp (lp ype) 7 8 = tan “(a Sao 
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R; = 
sing 
x 
2 
Ro = sind 
To a good approximation when 9 = 1/2, we have 
a= x, = Ri (1 - ~> =) (73) 
1 h_+h 
Sr 
Ry h 
b= R sin(Ad) Sin(Ad) 


The expression for <s*> is then introduced as a sSupple- 
ment to the direct path to the receiver in the following 
way: 


By utilizing the relationship 


re (74) 


for the intensity, the scattering intensity at the receiver 


and the intensity at the ensonified area are, respectively, 








7 8 a <s-> 
s ae ec 
nd < *> 
P, P1P1 
7 


Equation (70) can thus be written as 


AS(8) R? 


bd 
| 
- 


= 8m7R> <c}2> 
R2 
ea, = 282) 


2 
8TR2 <z'7> 
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Te 
mology T/T, = 10log 22 et 2 x 


ON: 2 


EB. = IL(R,) = K, (75) 


We then have to determine IL(R,) 


TL = SPL(1) - SPL(R,) = 20logR (76) 


1 


Semry) = IL(R,) = 10logt, 


SPL (1) SL 


ylelding 


IL(R,) = SL - 20logR, = K, (77) 


Then 


ie =K. + Ky (78) 


Suppl 
is the supplement to the direct path, and 


— sob 20logR, = K, 
See ee 
10 


Ihairect 


= anti log = K, 


ieee “rer 
ee 


iE = anti log a Ke 


suppl’ fret 
The total intensity of the receiver is thus 

ae it 
ii . direct + “suppl _, 4x =k 


i \ 5 
Tot’ ref L ef 


J. 6 


ILpo¢ = LOlogk, = Kz (79) 
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Thus, the effect of the randomly rough surface compared 


to the idealized free-field condition can be expressed as 


Be Tot direct ~ Ke ew 


A calculator program on a Texas Instrument 59 (later 
called TI 59) was developed to perform these calculations. 
A block diagram of the program is outlined in Fig. (26) 
and the programs steps together with a detailed description 
is given in Appendix B. The calculations are based on the 


following fixed data 


h =2m 
S 

h =o6om 
ag 

em 10 


Ww = 10 m/s (SS3) 
Then, varying the detection range from R = 2000 m to 100 m 


gives the difference between IL and IL Pilotaee qeaa 


Toe direct 
meg. (27). 


E. CONCLUSION AND DISCUSSION 

As seen from Fig. (27) the supplementary effect of 
the scattering from a rough surface in the high frequency 
case is negligible compared to the direct path. 

In saying so, it also should be pointed out that the 
Helmholtz-Kirchhoff's approach may be limited as it does 


not take into account the diffraction effects from crests 
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and throughs of the ensonified area, an effect which 
becomes increasingly important at low grazing angles, 
high frequency and when the rough surface is a super- 


position of swell and capillary waves. 
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VII. THE EFFECT OF SCATTERING AND ABSORPTION FROM THE SUB- 
SURFACE OCEAN LAYER 


A. GENERAL SCATTERING THEORY 

Generally when a region (volume) scatters sound, some 
of the energy carried by the incident wave is dispersed. 

The energy lost by the incident wave may be absorbed 
by the scatterers or it may be simply deflected from its 
original course. The amount of energy lost per second by 
the incident wave divided by the incident wave's intensity 
is called the total cross section o. of the region and is 


the sum of the absorption and scattering cross sections 





tt ale 
oo. = 
e if 
Pp 
where 
le = scattered power 
Ie = absorbed power 


The existence of gas bubbles in the subsurface ocean 
layer modifies the forward scattering in the following 
two major ways: 
1. The bubbles can resonate. When the bubbles are 
excited at a frequency near its natural frequency, 
it very efficiently absorbs and scatters ehesane2— 


dent wave. At resonance, the scattering and 
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absorption cross section of a typical bubble at sea 
1s of the order ee times its geometrical cross 
section. 
2. The bubbles change the effective COmpressibil lity Of 
the water and cause the speed of sound to be a 
function of frequency, i.e., the medium is dispersive. 
We will investigate and discuss these effects by separately 
estimating: 
1. The attenuation due to the bubbles and 


2. the refraction by bubbles. 


B. ABSORPTION MODEL 

The choice of model for the subsurface ocean layer depends 
on whether the medium has a teneous or a dense distribution 
of scatterers. When the bubble density is teneous, both 
"Single scattering" and "first order multiple scattering" 
approximation solutions are applicable. 

On the other extreme, when the bubble density is high, 
the so-called "diffusion" approximation can be used. 
Between these two extremes, multiple scattering effects are 
important. 

The multiple scattering theory, which in the limit also 
contains the first order approximation, will be used to 
estimate effects of attenuation due to bubbles, on the 
propagation of propeller noise from the target to the 


torpedo. 
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The geometry of the propagation model is illustrated 
in Fig. (28) where it is assumed that a plane wave is 
incident on a semi-infinite (disregarding the sea surface) 
slab of thickness x containing a number of randomly dis- 
tributed bubbles. The plane wave approximation is valid 
if the incident sound has a wavelength A much greater 
than the bubble-radius a 

ka<<l 
where 

kay /c=2T/d 
The receiver is located outside the slab and the beam 
pattern of the receiver is represented by the solid angle 
2. 

We are interested in the estimating of the total power 
received, taking into account the multiple scattering 
process in the inhomogeneous slab as well as the beam 
pattern of the receiver. 

The mathematical formulation of this problem is based 
on Twersky's theory of multiple scattering. Since the 
@meory is presented in Ref. 20 only, the basic formulation, 
major assumptions, and the end results will be presented 
here. 


The total intensity is the average of the square of 


the magnitude of the total field: 
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a e 
epee > = <l|<y™> + "| ?> = [<po> 2 + <|y."|?> (81) 


where 


the scalar field at the receiver location a 

see Fig. (29), is the sum of the incident wave 5 
and the contribution from all N scatterers. 
|<w*>|? is the coherent intensity based on the 
average field <p>. <|we>|?> is the incoherent 
intensity based on the fluctuating field Den: 

In Twersky's theory, the multiple scattering process is 
described by the following set of integral equations which 
Eq. (81) must satisfy: 

rele = levee [vee emeye|*>o(r_)dr_ 482) 


where 


a _ a a e -> 
ie = U. - UL Vie p(r,)dt 


is an operator representing all the scattering processes 

from s to a. (See Fig. (28).) It should be noted that 
Twersky's theory includes all the TUubeniig sea erenea waves 
that involve chains of succesSive scattering going through 
different scatterers. (See Fig. (29a).) However, the theory 
neglects the terms which include an individual 


scatterer more than once, as illustrated by Fig. (29b). 
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Thus Twersky's theory is expected to give good results 
when back scattering is insignificant compared to the 
scattering in other directions. 

As typical for most integral equations, Twersky being 
no exception, detailed complete solutions are difficult 
to obtain. However, Twersky gave an approximate solution 
merbea. (81) and Eq. (82), which according to Ref. 20 have 
been found to agree reasonably well with experimental data. 
This solution is based on the following main assumptions: 

-Backscattering is assumed to be small compared to 
scattering in other directions. 

-Scattering is mostly concentrated in the forward 
direction. This is reasonable based on the assumed 
random distribution of the bubbles; i.e., no rein- 
forcement of the radiation pattern occurs except in 
the direction of the incident wave. 

-The angle a is small, l.e., Oe = 0°, 

This leads to the following expression for the total 

intensity at the receiver: 

an |2> =T= exp (-90_x) exp (-90 _x) (83 


+q[l-exp (-po,x) ] 


where 
7 absorption cross section 
= scattering cross section 
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fl£l Pde 
g="r S 





ai lf] *aa. 
where 
f = the amplitude function 
eerie =<raction of total scattered power collected by 
eacereceiver as illustrated in Fig. 30 and 
© = is the bubble density, i.e., the number of 
scatterers per unit volume. 
For small values of ox we see from Eq. (82) that the 
coherent part dominates: 


mT = ~ (0, +0.) px (84) 


In this limit the multiple scattering result is equiva- 
lent to that obtained from single scattering considerations. 

For large values of px which corresponds to very dense 
or very wide slab of scatterers, the incoherent intensity 
dominates 


age = ing - J Ox Soy 


: ne > 
In this latter case, it is notable that when 2 =2T , 
the receiver collects almost all the scattered power 


ae (86) 
ree = Oo 0x 


The first case, representing the situation for teneous 
density of scatterers and/or narrow beam pattern of the 
receiver, gives a good approximation to the situation of 


interest in the thesis. It also represents the case for 
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which no scattered power is received. This will be 

approximately true for a narrow-beam width receiver. We 

will, therefore, Pee enigate the coherent intensity first. 
Also, the incoherent case, as represented by Eq. (86) 


will be investigated, where only losses due to absorption 


are incorporated. 


©. THE COHERENT INTENSITY CASE 
For the coherent case the intensity level after the 
incident wave has traversed a distance x 1S 


il 


gle exp[-(o,+0,) ex] (87) 


where 


i 
P 


The change in intensity over the distance x 1s 


incident plane wave intensity. 


i I. exp [-(0_+0_) 0x] 
AIL = L0log—= = 2 ____s$_= ___ (88) 
Pp i 


10logexp [- (9, +0.) 9x] 


The excess attenuation per unit distance due to bubbles 





is thus 
er Ama Ss ox loge 
: x 
= (o_+5,) 9 log e 
a = 4.34 0, 0, in dB/m. oe 


However, this only takes into account bubbles of one 
size. In a bubbly medium there is a spectrum OL adil. 


fie 





The probability density function for ELNGLIGea SUbple Size 


between radii a and a+da is 


W(a) = — (90) 


where 


Cc 


f W(a)da = 1 
O 


n(a)da is the number of bubbles per unit volume 
having radii between a and a+da. It is common to 


use da=l1 um. 


0 = n(a)da is the total number of bubbles per 
peaks volume. 
As the extinction cross section also is a function of 
the radius (See Eq. (93).), the absorption due to 
bubbles is obtained by integrating Eq. (89) over all 


possible radii 


co 


aw 34 i) a, (a)n(a)da (91) 
O 


To calculate the absorption coefficient, the extinction 
cross section o, must be derived from the general bubble 
dynamic relationship. This is done in detail by C. Clay 
and H. Medwin [Ref. 21] from which the following results 


are taken. 
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The scattering cross section 





I 
= oe) 
P 
= 4na* 
Gee =) =| +6" 
where 
2YbgpP 
~ eee A, *3 
£. resonance frequency = Se o, ) 
f = fo = operating frequency 
6 = damping constant = Or + by + 6, 
£ 
= Gy ene 4u 
= ka + (5) (=) + eae 
A 
a = bubble radii 
¢ _ 3 (y-1) (20822 h x + Sin A) -2(cos h X - cos x) = 
X*(cos h X - cos X) + 3(y-1) (sin h X - sin X 
Sy ee 
i a (—P2 25) % 
i: g 
= = thermal conductivity of gas 
Za 
= — — ~ Z 
oe density of gas Pgallt (yay! OzZ) 
eee o Ley of gas at sea level 
tT = surface tension 
Pp, = 1.013 x 10°(1+0.1z) 
Za oubpDle depth in m. 
Cae specific heat of constant pressure of gas 
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5 
ll 


Shear viscosity of water 


7/5, for diatomic gas 


=a 
il 


oS 
i 


A density for sea water 


= 


b = [1+ (S) 217? [1 = =A ie Sin h pe = sin X 
X cos h X - cos X 
2T lL 
oe et ——— (1 - ——) 
Pia 3yb 
Furthermore: 
2 A 
5, = J,t0, = __4ma* (6/ka) (93) 
Ga 2): -1]7+6? 
and 
J. = C 708 (94) 


A detalled computer program, as outlined in Appendix B, 
was developed for the TI 59 to handle the derivation of 
Jos Oy! and o based on an assumed receiver depth of z = 6 m. 

For the 60 kHz case, both the extinction cross section 
J. and the absorption cross section Jd, are Given iiebaa. 3 o 
as a function of bubble radius a. 

Similarly, Fig. 32 gives J. and J. for the 30 kHz case. 
Superimposed on these figures are the curves for n(a)da as 
calculated from the following: Figures 5 - 7 of Ref. 2 
give the resonant bubble densities in a1 um band as a 


function of depth and with the wind speed as parameter for 


the three discrete frequencies 12 kHz, 38 kHz and 120 kHz. 
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Based on these data, Fig. 33 shows the interpolated bubble 
density as a function of resonance frequency for sea state 
m™ 3S, and 6. | 

Furthermore, both A. Lovik [Ref. 2] and H. Medwin [Ref. 3] 
found that the bubble density function n(a) decreases with 
increasing bubble radii as 

naa * 
where H. Medwin [Ref. 3] found the power law of: 

xX = 4 for a<50-80 um. 

x = 2 for a>50-80 um. 
and A. Lovik [Ref. 2] found the power law of: 

x = 4.2 between 38 kHz and 120 kHz 

x = 2.6 between 12 kHz and 38 kHz 
which, averaged over the depth interval, corresponds to the 
bubble radii of 380 uwm(12 kHz), 120 um (38 kHz) and 49 um 
(120 kHz) . 

As suggested by A. Lovik [Ref. 2], the discrepancy 
between the two obServations is not great and may be due 
to the few measuring frequencies used in the work of 
ma. LoOvik . 

In summary, Fig. 33, from which we obtain the appro- 
priate resonant bubble density ina 1 um band n (ap) 


together with the power law naa, comprise the full 


knowledge of n(a)da. 
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Performing a multiplication of oe and n(a)da, we 
obtain Fig. 34 and Fig. 35 for the 60 kHZ and the 30 kHz 
case, respectively. 


Pamally, the integral 


oO 


/ J, (a)n(a)da 


was evaluated using a numerical integration based on 
Simpson's discrete approximation programmed for the TI 59 
and documented in Appendix C. 

Based on the above, the following absorption coeffi- 


cient for the coherent case is obtained for f = 60 kHz 


oO 


a = 4,34 i J, fa)n(a)da = 4.34 (2.016 x 107") 
O 
7 ib 
=e5./5 xX 10 dB/m. 
a * 0.88 d4B/m ee) 


with 
n(ap,) = 1000 
n(a) a a-4 
Z = 6 m@. 

Thus, it is seen that the attenuation due to bubbles 
is considerably greater than the "normal attenuation" 
due to chemical and viscous relaxation processes in sea 
water, which for the 60 kHz case is 


fo On02eaB/m. 
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Thus, the total absorption coefficient for the 60 KHZ 


case is 
on 0250 dBm. (96) 


For £ = 30 kHz, the absorption coefficient due to bubbles 


in the coherent case is 


oe) 


a = 4.34 yf aa) ital) da r= 43 4a3 188 75m x 107°) 

0 

= 1.665 x 10! aByn. 
a = 0.17 dB/m (97) 
with 
n(a,) = 20 
n(a) a —— 
Z= 6m 
Zz. 


For this case, the power law dependence of n(a) aa 
gives 
a = 1.655 x 1077 
Thus, the difference in power law dependence makes no 
Significant difference in the absorption coefficient. 
The absorption coefficient due to chemical and viscous 
relaxation processes is at 30 kHz 
e =80.012 dB/m. 
The total absorption coefficient in the coherent 


case is 


a = 0.18 dB/m. (98) 
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See tit, LENCOHERENT INTENSITY CASE 

For the incoherent case, where only losses due to ab- 
sorption are included, the product om and n(a)da for the 
60 kHz and 30 kHz cases are given in Fig. 36 and Fig. 37, 
respectively. Performing a numerical integration based on 
the Simpson's discrete approximation leads to the follow- 
ing results: 


For £ = 60 kHz, the absorption coefficient due to 


bubbles is 
7 7 -1 
a= 4.34 f a (a)n(a)da = 4.34 (1.66 x 107°) 
0 
a = 0.72 dB/m. (99) 


with 


n (ap) = 1000 


mta) a a 
z= 6m. 
Adding the "normal attenuation" in sea water for 


f = kHz, yields a total absorption coefficient of 


a = 0.74 dB/m (100) 
For £ = 30 kHz, the absorption coefficient due to 
bubbles 1s 
a = 4.34 /0 (a)n(a)da = 4.34 (2.975 x Oma CO) 
0 
a ~ 0.13 dB/m 
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Adding the "normal attenuation" in sea water for f = 30 kHz 
ylelds a total absorption coefficient in the incoherent case 
OL 


a + 0.14 dB/m fale) 


E. SUMMARY AND DISCUSSION OF THE BUBBLE ATTENUATION 

As seen from the above, the attenuation due to bubbles in 
the subsurface ocean layer 1S important for high frequency 
and high Sea State. 

Accounting for the effect of bubbles at Sea State 3 in 


summary we found the following absorption coefficients in 


adB/m: 
f = 60 kHz fo — eo Oe s 
coherent incoherent coherent incoherent 
a = 0.88 on = Oye eee (e243 


The main assumptions were: 


-the back scattering is small compared to that in 


other directions. 


-the scattering is mostly concentrated in the 


fomwarad direction. 
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-the angle oe between the scatterer and the 
receiver is small, i.e., era 40 oe 
-the depth of the receiver is z = 6 nm. 
-the sea state is 3. 
The range dependent portion of the passive sonar 
Eq. (1) TL=-20logR-aR for both the 60 kHz and 30 kHz 
scattering results can now be compared with the reference 
@ata obtained in Section V. 
For £ = 60 kHz, Fig. 22 gives TL=-20logR-aR as a 
meumctioen of R with 
ae © =—sOnse71aB/m, the “normal attenuation" due to 
chemical and viscous relaxation processes. 
2. a= 0.9 dB/m, the total attenuation including the 
effect of bubbles in the coherent case. 
3. a= 0.74daB/m, the total attenuation including the 
effect of bubbles in the incoherent case. 
For the same example as in the reference model, Fig. 22 


yields the detection ranges for 60 kHz. 


R = 1200 m for a = 0.021 dB/m | (103) 
ne oemm for «a = 0.90 dB/m 
R = 70 m for a = 0.74 dB/m 


Not surprisingly, this result seems to exclude the possi- 
bility or having both a searching and attack depth near 
the surface, i.e., z = 6 m, for a torpedo system operating 


in a high frequency region, f£ = 60 kHz. 
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Similarly, for f = 30 kHz, Fig. 23 gives TL=-20logR-aR 
ae a function of R with 

1. a = 0.012 dB/m, the "normal attenuation" in sea 
water. 

2. &% = 0.18 dB/m, the total attenuation including the 
effect of bubbles in the coherent case. 

3. a = 0.14 dB/m, the total attenuation including the 
effect of bubbles in the incoherent case. 

For the same detection example as in the reference model, 


Fig. 23 yields the detection ranges for 30 kHz: 


R = 2400 m for a = 0.012 dB/m (104) 
R= 250 m for a= 0.18 dGB/m 
Pee ot0O m for a= 0.14 dB/m. 


Again, the bubbles give a major decrease in the detection 
range. A detection range of R = 250 m seems marginally 
acceptable as the turn rate requirement for a pursuit 
homing trajectory may become excessive. 

The above results are summarized in Tables TIIand IV 


for the 60 kHz and 30 kHz, respectively. 


P. THE REFRACTION BY BUBBLES 

The presence of bubbles in the sea water affects the 
speed of sound (phase speed) primarily because of the change 
in compressibility. The derivation of this dispersive rela- 


tionship on the sound speed has been done by H. Medwin 
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[Ref. 22]. H. Medwin showed that the bubbles with resonant 
frequencies greater than the incident frequency decrease 
the sound speed, while bubbles with resonant frequency lower 
than the incident wave increase the sound speed. 

Furthermore, H. Medwin [Ref. 22] predicts the sound speed 
gradient due to bubbles as a function of depth for a frequency 
range and wind speed compatible with our domain of interest. 


He found the gradients 
i 


g = 3c/3z = 0.26 Ss at z=On. 
g = dc/faz = 0.016 st at z = 10 nm. 
g = dc/az = 0.005 st at z = 20 m. 


For comparison, the sound speed gradient due to pressure 

in an isothermal water is 
g = 0.017 st 

This shows that the rays in the top 10 m are strongly 
influenced by bubbles. However, with respect to our surface 
Scattering model (Section VI) where both the source and the 
receiver are situated very close to the surface and where 
relatively short propagation distances are encountered, 


this refraction phenomenon is assumed to have negligible 


effect. 
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ViTI. THE TURN RATE LIMITATION 


AS pointed out in the previous section, the presence 
of bubbles near the surface may significantly reduce the 
range at which the target can be detected. 

In this section, the turn rate necessary during 
pursuit homing at the previous estimated detection ranges 
will be studied. A computational procedure will be used 
to determine the range of angles on the bow (AOB) of the 
target at the beginning of homing which lead to miss on 
the initial attack. A trajectory where the torpedo velocity 
vector always is directed towards the instantaneous target 
position is called a pursuit homing trajectory. The 
derivation of the pursuit homing trajectory follows 
P. van Nostrand [Ref. 23] and is based on the geometry of 
Fig. 38, where 

r = approach angle, i.e., angle between ship 
velocity vector and the line of sight. 
AOB = 180 - $6 

The equation of motion is obtained by taking the 

component along r and the normal to r, yielding 


— ~ 105 
= ee cos $¢ Vin ( ) 


rd = -v. sin ¢ (106) 
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where 


r= range rate 
¢? = turn rate 
We = target speed 


Vin torpedo speed. 
Prot aing eq. (105) by Eq. (106) yields: 


et 


Te aie EOEyO) Oe AO 


KY [rie 


and defining 


Yr A 
V = —p 
S 
yields 
ee Po 
=F Pad cot $)9 (107) 
P-1 
= eisai ) 
aa ‘(1 +c0s86) a 
ian, ancroducing the initial conditions: toe 6 where 
ro? initial detection range 
oes initial approaching angle 
yields p 
(l+cos Oe) 


; p= 
aan ed 
For the geometry of Fig. 38, the turning rate is given by 


Bq. (105). 
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Vy 

$ = --=> sin ¢ (110) 
Pu@eettucing Fao. (108) into Eq. (110) yields: 

; V,(ltcos o)® 

ee SS ioa.6°° Cl) 
From Eq. (108) we see that r~0 as 970, i.e., 6 tends to 
zero as the torpedo approaches the target ship. It is 
further of interest to determine the limiting value of the 
turning rate as the torpedo approaches the target for 


various values of the parameter p. This is done by taking 


the derivative of 6 with respect to ». Thus, from Eq. (109) 


we get 
ai Vs 1-P P 
_ = - (sin 9) [lt+cos¢]” [2 cos ¢ - p] (112) 


Bor l<p<2, we see that Eq. (112) is zero at 
cos 6 = P/2 chigs 
pe es - (2/2) 


with an associated maxima 


ey 2 


V_({1+P/2)" [1-(P/2)7] 


Omax K (i 








Furthermore, the turn rate for p<2 at impact 1s zero as 


eeeeriat of Eq. (109) yields: 


V 


lim $= - — lim(sin ifn tdteos o) So Gini) 


o+0 +0 


oi 





For p=2 we see from Dq. (111) that 


V 4V 


lim $= - 2 lim(1+ cos Dio ee 2 (116) 


p70 g+0 
For p>2 we see from Eq. (111) that there is no zero 
malue of do/dae between 6 = 0° anda o = ‘Leven Since all terms 
eemeag. (112) are nonzero terms for Or cane0- Furthermore, 
the turn rate for p>2 at impact is 
Yr ..  (1+cos o)? 


aaa oo, = |= —— Lam = 


53 | +00 (ete) 


ono (sin 6) 


as (Sin ?P°Since (p-2)>0. 


Furthermore, as 6+180°, we get from Eq. (112) 


V 


lim =~ = =e (2+p) aim (sin On (icos )? 
97180 
Evaluating 
ieee ailiscos 6) =edim (sin 0) (=EgeS 2)” 
o+180 dae 0 


by applying 1' Hopitale rule to the term 


Weapligeos ss — >./=5am > 
lim 35h =e lim Soom 0 
6+180 6-180 


1+cos ONE 


approaches zero for 
sin 


Thus, the product sin 9{ 
any p>l, since both terms in this product approach zero 


as 9-180. 
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Figure 39 shows a plot of computed values of |bK/V_ | 
for different values of p>1. From this we can draw the 
following summarizing conclusions. 

bem !<p<2 the turn rate has 

-a maximum value at $=cos7)(p/2) 
-the zero value at ¢=0° and =180°. 
Benuep=—2 the turn rate 
-is zero at o=180°. 
“monotonically increasing with decreasing 9 
approaching the value (4V _/K) when >=0°. 
For p 2 the turn rate 
-1s zero at >=180°. 
-is monotonically increasing with decreasing 6 
approaching © as o=0°. 

As seen from the above for p>2, the turn rate increases 
monotonically with decreasing approaching an infinite 
PUeaetate tO hit a point target. To avoid this singularly 
we must make some provision. If the torpedo's maximum 
turn rate is exceeded only at some very small range, a hit 
is likely. 

Figure 40 illustrates the hit criterion used. Assume 


a rudder of length L, is situated directly behind the 


R 

propeller which is idealizer as an acoustic point source. 
If the torpedo becomes turn rate limited at some range 

r=r', it at best can proceed along a circular path which 


lags the desired pursuit trajectory, or at worst it can 


= 





Me@ee acoustic contact and go into "hold-in," maintaining 
a constant heading at the angle ¢'. In this last case, it 
will cross the line of target ship's course in the time: 

= ae (118) 
If this advance is less than the length of the rudder Lp: 
the torpedo will impact the rudder. Hence, the limiting 
Genedttion is 

gt 


p ~ OR 


ae pL, (aS) 


For this analysis, we have arbitrarily chosen 


oS = 03 on (12 0) 


Thus, the torpedo's turn rate has not been exceeded 
when the range to the target is 


r' = pL, = 3 Pp (eZ) 


and a hit is assumed. 

Now we can analyze the cases p>2 and 1l<p<2 on the same 
footing. If we match the torpedo's maximum turn rate to 
a particular range r', it is certain that Hee Sane one isle 
has not been exceeded earlier in the run. 

This turn rate-range matching is done by substituting 
the "hit-range" defined as r' = 3 p into Eq. (106) and 


Salgmmag it for sin’ >, yielding 
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V 


6 =e S : 
max 3p Sin d 
b. n( 3p) 
: 2a 
Sino=- 
= 


Since 9? is always negative, sin 6 is always positive and 
equal to 


3p| | 


ee max 


V 
S 


(le ) 


Two values of $ satisfies this equation, and they are denoted 


2 3p) ° 
¢, = Sin I oe for ¢< 90° i223) 
S 
= OS 
d, = 180 o, 


A unique value of K may now be found using Eq. (111) for 


each of the angles On and pn: 


V.(l1+cos ae 
7 





| >| (sin 0)? 


Max 


These values of K are designated as Ky and Kay and for 
each there is a corresponding value of ie Prom cg 0 (109) 


with ® given by the initial detection range 


eos and (Op 


In the case of l<p<2, the turn rate does not necessarily 
increase monotonically during the pursuit homing trajectory, 
and we must check that either On or Ons respectively, is not 
beyond the angle corresponding to Cane 


oe 





As seen from Eq. (113), the turn rate reaches a maximum 


at an angle given by 
pee = cos 1 (p/2) 


with a corresponding turn rate given by 





| V_(1+P/2)P [1-(p/2)2y17P/2 
ale 7 K 


digas, the value of K = K* for which the limiting turn 
rate is achieved 


_ (L+P/2)* [1-(p/2) 2] °7P/? 


K = - 
Pl maxx 


iwtoerimportant to note that, if 9 Or) the "A" 


a! 
trajectory is invalid since the maximum turn rate of the 
torpedo, reached at the range r = 3p at 9 a 1s exceeded 
earlier in the trajectory. In that case, the limiting 
trajectory is the "*" trajectory. Along the same lines, 
we argue that the "B" trajectory always is a limiting 
trajectory, since the equality rs (P dp cannot be 
satisfied. This follows from the fact that 

O°<o* < 90° 


and Pe 90° as ? is the supplement of Pare 


Below some critical ship speed, the torpedo will not 

be turn rate limited. This speed is obtained when the 
; . O 
trajectory is normal to the ships velocity vector, 9=90 , 


at a range of r = 3p. 


96 





Thus, from Eq. (106), we obtain 


: _ O 
(3p) 1o1 as Nee by) FN 
yielding 
(Vour = (3P) lol, (124) 
Then introducing p = Vid (Ve) ae we get 
toe. = [(3ve1o| }3 (125) 
Ss NL aL max 
Furthermore, we see from Eq. (108) that for >=90° 
(K) uy = 3P (126) 


Again, by probing Eq. (88) we can obtain the corres- 
ponding values of een designated (9 ane where the subscript 
"NL" is used to indicate the "no-limit" boundary point. 

The computational procedure is based on: 

1. Assuming a torpedo speed V = 35 kts. 

2. Using an initial detection range (beginning of 
homing) ES 200m. 

3. Assuming the following maximum turn rates: 
ave. siatey Ac! ba sz o7 Ss C2 16. (67s 
Cl 2 eyes e. 36 0/s f. 48 o/s 

4. The following range of target speeds 


0<V<25 kts 


The aim of the computation is to determine whether the 
limiting Pit and hence AOB, is governed by $* or by the 


turn rate at the range r = 3p. 
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The calculations are devided into two parts, and are 


Bexrtermed om a TI 59 calculator. 


For a given ship speed, Part I gives the sequential 


O., K 


Galculiations of Par K O* and K* for each of the 


Pua By ae 
turn rates. The computer program is given in Appendix E. 
Then, in Part II, the probe calculation for (ds) ar 
(90) 5: On)” = (o5)a Be are performed. The program is 
given in Appendix E. 
The limiting results are given on a polar plot, Fig. 41 


wath 


— O _ 
(A0B), p,* — Eee (05) a,B,* 


As seen from Fia. 41, for a target speed of 15 kts, we 
need AOB>85° at 6 for a maximum turn rate 16 o/s in order 


to have a hit at the first attack. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


The detection performance of a passive homing torpedo 
used against shallow-draft surface ships operating in 
moderate sea states was investigated. Attention was 
focused on the effects of scattering from the randomly 
rough sea surface and scattering and absorption from the 
bubble-dominated inhomogeneous layer just below the sea 
surface. The effects of these two scattering mechanisms 
were separately estimated and their relative importance 
were compared. 

The passive sonar equation was used to predict the 
performance of the homing system, and the detection range 
considering these two scattering effects was obtained and 
compared to the detection range based on a reference 
model. An idealized propagation model was used as 
reference of comparison. This reference model was based 
on a noise source model for a cavitating propeller, the 
operational characteristics for a square-law detector, 
and a transmission model associated with a homogeneous, 
unbounded medium. 

Due to high frequency, moderate sea state and low 
grazing angles, the scattering from the randomly rough 


sea surface was found to be dominated by the direct path. 
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Maes result includes effects from geometrical 
shadowing. 

The effect of scattering and absorption from the 
bubble~dominated, inhomogeneous subsurface layer was in- 
vestigated using multiple scattering theory. Both the 
coherent and incoherent limits were investigated by incor- 
porating the associated absorption coefficient into the 
transmission equation. 

The effect on the sound speed from the bubble content 
was found to be negligible. 

At the assumed depth setting of 6 m for the torpedo's 
search and attack phase, the scattering from the bubbles 
increased the transmission loss. This increase depended 
on the frequency and the wind speed. 

Two operating frequencies were investigated, 60 kHz 
and 30 kHz. For both cases, bubbles significantly 
decreased the detection range. 

For a torpedo system operating at the high frequencies, 
e.g., 60 kHz, the result indicates the inadvisability of 
using a searching and attack depth near the surface, l.e., 
Z>= 6M. 

For an operating frequency around 30 kHz, the calculated 
detection ranges is such that the turn rate requirements 
Een a pursuit homing trajectory become excessive. For a 


maximum turn rate of 16 o/s, this limitation can be avoided 
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by adapting a tactical procedure where the angle on bow 
(AOB) at the beginning of the torpedo attack is greater 
than 85°. 

At sea state 3, the results show a consistent and 
general trend towards the need for lower operating fre- 
quency in order to increase the detection range. An 
operating frequency below 30 kHz seems indicated. 

Furthermore, a search depth below the bubble-~-dominated 
subsurface layer, i.e., z>15~20 m would result in an 
increased detection range. 

To reduce the operational limitations induced by the 
scattering and absorption effects, a high maximum turn 
rate together with a variable speed capability, where 
p=2 would be beneficial. 

The result of this analysis has clearly demonstrated 
the importance of environmental factors on the torpedo 
capability, and is useful in giving insight into the be- 
havior of a homing torpedo during its critical attack 
phase. 

A valuable follow-on of this study would be an investi-~ 
gation of the effects of the bubble~dominated subsurface 
layer on target validation and pitch plane steering when 


the torpedo search depth is 50-60 m. 
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APPENDIX A 


DETAILED OCEANOGRAPHIC BACKGROUND MATERIAL 


A. GEOGRAPHIC DESCRIPTION 

The Norwegian coastal waters constitute the eastern 
boundary of the Norwegian Sea. Although some general 
aspects related to the Norwegian Sea will be covered, 
this analysis will be concentrated on the Norwegian 


coastal waters above 68°N. 


SB. i WwitND 

The northern region of the Norwegian Sea is affected 
by the Polar Easterlies and the southern region by the 
prevailing Westerlies. There are two dominant air masses 
which are relatively permanent: 

-The Greenland high. 

-The Iceland low. 
These pressure systems produce storms which are carried 
across the Norwegian Sea ina belt from Iceland towards 
the Norwegian Coast causing steady precipitation and wind 
most of the vear. The steep Norwegian Coast has a con- 
Siderable influence on the winds and consequently also on 
the waves in the coastal waters. The main general modi- 


fications are that the streamlines tend to run parallel 
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to the coast and that wind and sea increases with distance 
from the sheltered coast into open ocean. 

Strong local variation may occur. The most important 
of these are the marked local increase in wind speed in 
areas where the coast sharply changes direction. One such 
"corner effect," caused by the confluence of the streamlines, 
occurs near "Nordkapp" (North Cape). Also of importance 
are monsoonal effects due to the different heat capacities 
between the ocean and the continent. Drainage of cold air 
from the inland valleys in the wintertime causes a marked 
increase in the wind speed in several areas along the 
coast. Most of these coastal effects are significantly 
dissipated at distances of approximately 50 nmi. from the 
coast. The fact that the wind tends to blow along the 
coast is clearly demonstrated in Figs. 12 and 13 which 
include data from weather stations from "Hillesoy" to 
"Ona" and "Myken" to "Furuholmen" respectively. The high 
frequency of offshore winds is caused by the drainage of 
cold air from the inland valleys during winter time. 

That this phenomenon is closely connected to the coast 

is illustrated by the fact that it is missing at weather 
Station "Skomvaer" situated approximately 50 nmi. off the 
main coast. A frequency distribution of observed wind 
speeds along the coast, obtained from Ref. l, is presented 


in Table V. A summary of this table follows: 
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TABLE V 


PREQUENCY DISTRIBUTION OF WIND SPEED IN PERCENT PER YEAR 
AT WEATHER STATIONS ALONG THE NORWEGIAN COAST 


m/sec 0-7 8-13 14-20 =21] 
Beaufort 0-4 5+6 7+8 29 
Ferder 67 28 5 OL 
Lyngor 79 19 2 Oe2 
Lista Sr 25 5 O22 
Utsira 13 22 5 02s 
Hellisoy aa 1S 4 OiZ 
Krakenes SiS 28 eA 2.7 Jan-Dec 
Ona 70 22 7 ae ee 
Sula 61 Be ical 
Nordoyan 49 a5 14 Ie 
Myken 64 26 9 i lehag” 
Skomver 58 oil sO O29 
Andenes 712 18 Ore 
Torsvar Ue 23 0.4 
Fruholmen 54 Be 12 ies 
Vardo iS ZZ O.1 
Biornoya 63 BA Onaga > = 197 5 
Polarfront 46 30 14 1.4 
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The highest winds are reported from the areas between 
62°N and 68°N. The trequency Of@hiohuwinds inerhis area 
are significantly greater than those found at stations to 
the south of 62°N and also to the north of 68°N. The wind 
condition around North Cape, are very severe. This area 
can be compared with the other coastal area of high wind 
speed such as "Stadt" (represented by the weather station, 
“Krakenes"). In these areas the frequency of storms is 
greater than at "Polarfront" situated in the open Norwegian 


Sea. 


G WAVES 

Frequency distribution of significant wave heights are 
represented in Table VI. 

The station north of 68°N is characterized by compara- 
tively small frequency of high waves. Even at "Furuholmen" 
where wind conditions are very severe, the frequency dis- 
tribution of significant wave heights is similar to more 
sheltered areas like "Utsira." 

The seasonal variation for the area of interest around 
70°N is given in Figs. 3 and 4. The average monthly dis- 
tribution of significant wave heights for a typical 
station like "Andenes" is given in Fig. 3 and Fig. 4 


gives the yearly distribution for this station. 
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TABLE VI 


FREQUENCY DISTRIBUTION OF SIGNIFICANT WAVE 
HEIGHT IN PERCENT AT THE WEATHER STATION "ANDENES" 


STATION: Andenes YEARS: 1949-1972 
SS 0+1+2 3 4 5 6 7 8 9 
Ho 0 Seen Demeel. Ss — 245-4250 ~6.0) = 9.0—. 13.9 N 

0.4 Lee 2.4 a9 «Seo: UB ee)6~6l(lUd1L3.S 
JAN 15.3 Zon es oe 15S 0.18 0.00 2542 
FEB ey] woe bl. 8 lhl 0.14 0.00 2295 
MAR 18.1 28.9 30.8 14.3 5.8 : OF? 6 OF0G™ 2529 
APR 2a 5 Boe be 27245510.8.. 3.4 Peo O4 2293 
MAY 30.9 Boe5 23-9 ee oO On 1720 
JUN 37.0 36.9 20.2 4.8 0.00 L565 
JUL 42.6 wiees | 6LL 7. Bo7 ; 0.00 1618 
AUG 40.6 36.4 18.0 Ceo a9 eee 1621 
SEP 29.7 Bae 3 ecto 8.9 6701 2301 
OeT 21.8 eis 62954 «=o 24 ; OF; OL 2o74 
NOV oar 7 Be Seca | 6L le 2 1.1 0.06 2303 
DEC o> Berl 86h3 22°C Oey C2 C008 2379 
YEAR 26.6 Be 5 wee Oo Ole. 3 - Oo 0606) -0.00 

N = Number of Observations 

SS = Class Interval, State of Sea 

Ho = Significant Wave Height 


106 





D. AMBIENT NOISE 

Due to the fact that the Norwegian Sea is physically 
separated from the Atlantic by the Faerdoy-Shetland-Iceland 
ridge, and from the Greenland Sea by the Jan Mayan ridge, 
little long distance shipping noise is transferred into 
the area. This, combined with relatively low shipping 
traffic in the central and northern parts of the Norwegian 
sea, produces a relatively low ambient noise level for the 
frequency band 100<f<1000 Hz, especially when noise from 
marine life is not included. Further, in the absence of 
nearby shipping and marine life, the ambient noise level 
in the frequency band 1<f<50 kHz is, according to Ref. 4, 
dominated by the wind. 

A typical area for the central part of the Norwegian 
Sea can be represented by the weather station "Polarfront" 
at 66°N, 2°E. Reference 1 shows that there is approximately 
a 15% chance of finding wind forces of Beaufort>6. These 
effects predict a moderate ambient noise level in the 
frequency range 100 Hz-50 kHz for the central and northern 
regions of the Norwegian Sea. 

The shallow coastal Norwegian waters are, according to 
Ref. 4, typically 5-10 dB noisier than the corresponding 
deep waters. However, great variability caused by local 
traffic, fishing fleet activity, marine life, and local 


wind conditions makes ambient noise level prediction 
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difficult in these areas. This means that accurate ambient 
noise level determinations have to be made on the spot, as 
it is both site and time dependent. 

A noticeable influence on the ambient noise level is 
rain, which is a year around feature along the Norwegian 
coast. As seen from Fig. 8, taken from Ref. 4, rain has 
a tendency to produce a constant high ambient noise level 
over a large frequency range, thus dominating other 
effects. Furthermore, for the upper frequency of interest, 
1.e., around 60 kHz the lower bound for the ambient 
noise is determined by the thermal agitation, see Fig. 8. 

In determining the figure of merit (FOM) for a 
passive sonar system, the noise level will be the larger 
of either the self noise or the ambient noise. Fora 


torpedo the self noise will typically be dominant. 


=. SOUND SPEED PROFILES 

Again, concentrate on data relevant to Norwegian 
coastal waters. According to Ref. 6, which covers the 
southern part of the Norwegian coast, low sound speeds 
are common because of the influence of water from the 
Baltic Sea combined with fresh water drainage from the 
fjords. Furthermore, great variability, both seasonal 
and within seasons, is encountered. Figure 9 obtained 


from Ref. 5 gives a picture of the sound speed profiles 
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Bor the northern Norwegian coast. Again, large variations 
are common. Noticeable in both sets of data is a typical 
seasonal pattern of strong cooling of the surface layer 
during winter and a similarly strong heating during 
summer. Furthermore, note that the minimum and maximum 
are relatively shallow, i.e., less than 50 m. 

Also characteristic is the influence of the cold and 
fresh melt waters drained out through the fjord-arms 
during spring and summer. 

To illustrate the sonar problems associated with these 
sound speed profiles, ray paths for the extremes of Fig. 9 
are shown in Figs. 10 and 11, where the source is 3 m below 


the sea surface. 
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APPENDIX B 


SURFACE SCATTERING TI 59 PROGRAM 


& INTRODUCTION 

This program gives specular scattered power at the 
receiver versus incoming power at the randomly rough 
surface in the high frequency limit according to Eq. (70). 

Shadowing of surface areas by other parts of the 
boundary are taken into account by the bistatic shadowing 
Pamecion S$(6). 

Furthermore, the program gives the effect of the 
randomly rough surface compared to the idealized free- 
field condition as expressed in Eq. (80). The results 


of these calculations is given in Fig. 27. 


B. PROGRAM STEPS 

A block-diagram of the computer program is given in 
rig. 26. 

The program uses the partitioning ratio of program 
moemcdeta spaee according to code 4 OP17. The users in- 


structions are as follows: 
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Procedure Enter Press Display 


Enter data Detection range 2nd A a 
Enter data Source depth RS Be 
Enter data Receiver depth R/S ee 
Enter data Beam width R/S oe 
Enter data Wind speed R/S wt 
Calculate 4, ae 2nd B A 1 
BR, AR2,A (6,Ri,R2,A) 
Gateutate <z'*> -- 2nd Cc <ri2>t 
Cofkeullate v == 2nd D vi 
Calculate erfc v -— A erfc vi 
Calculate S(96) —— B s(e)t 

a+ > ég2>t 
eamcarate <p,p,*> ee G gee oe 
Calculate AIL ae D aru 
fe 


Meese values are printed automatically if the calculator 
is connected to the PC-1LOOA Print Cradle. 


For the error function complement we have 


er 
erfc(v) = A e * da 
V 


in contrast to the normal distribution 


ats 


>) meee 
my = == f{ ae /* at 
Y27 a 


However, there is a linear relationship between the 


ie functions. 
lel 





is a modified program from Texas Instruments 


The numerical equation used to calculate the erfc(v) 


[Ref. 


= Z(v) [biat+b2a*+b3a%+bya*+bsa°] 


erfc(v) 
where 
eS = .231649 
be = 451673691 
be = -.504257336 


Pepe 2 51939026 


Ep —2. 565934662 3 


Ds 


eo 292139 


24]. 


The program steps are listed below, giving location 


(LOC), code 


000 
001 
OZ 
003 
004 
005 
006 
007 
008 
009 
010 
O11 
O22 
013 
014 
02> 
016 
017 
Ors 
OdeS 
020 
O21 
022 


7o LBL 
16 A’ 
ae ST 
00 00 
2 PRT 
ol R/S 
a2 StO 
Ol Ol 
oy PRT 
onl 3R/S 
42 STO 
Oz 02 
ge PRT 
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#2 STO 
03 03 
Sore Rt 
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42 STO 
04 04 
Ge" PRT 
Sr R/S 
76 LBL 
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O27 
028 
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O31 
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O25 
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OS7 
038 
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043 
044 
045 


(COD), key symbol 


ipeZ 


(KEY), 


B! 

( 
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RCL 
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Bee ENDix Cc 


BUBBLE DYNAMICS TI 59 PROGRAM 


A. INTRODUCTION 

This program gives the resonant frequency f. and the 
damping constant 46 for bubbles according to Eq. (92). 
Furthermore, the program gives scattering cross section 
Oo., extinction cross section Oa: and absorption cross 
Section o, asa function of bubble radius a, incoming fre- 


quency £, and depth z. The results of the calculations 


are given in Figs. 31 and 32. 


B. PROGRAM STEPS 


The user's instructions are as follows: 


Procedure Enter Press Display 
Enter data Bubble Radius 2nd A at 
Enter data Incoming Frequency R/S gl 
Enter data Depth RS zt 
Calculate xX -- 2nd B xt 
Calculate l 
coshX and coshxX 
sinhx -- Zine «C sinhxl 
Calculate d/b -- 2nd D a/bt 
Gamculate b -- 2nG & bt 
1 
Calculate 8 id A 8 
Calculate f == B f l 
r r 
iE 
Calculate 4 -- Cc § 
Calculate o -- D oe 
Calculate - 1 
and o -- E oe 
a 
cee 
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1 
mese VYedues are printed automatically if the calculator 


1s connected to PC-lLOOA Print Cradle. 


The program steps are listed below giving location (LOC), 


code (COD), key (KEY) and comments. 


000 76 LBL 045 55 = eG O01 On 
001 16 A' 044 43 RCL 087 99 PRT 
G92 82 STO O45 i il 088 91 R/S 
es 11 11 046 54 ) 089 76 LBL 
gee 99 PRT 047 85 + 090 18 Cc! 
te5 91 R/S Os Fe.4 O91, 753 4 
Geo 42 STO 049 00 0 092 53 ( 
oo7 b2 12 050° 00 @ 093 43 RCL 
008 99 PRT Csi Of 1 094 01 01 
009 91 R/S O52 02 2 O95 22) INV 
mo 82 STO 053 09 9 096 23 LNX 
Ge 86s «sd 3 054 65 x 097 85 + 
012 99 PRT 055 53 ( 098 43 RCL 
fs «20691 CUR/S 056 O11 099 O01 Ql 
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Gms 1/7 B' 058 93 - TO 22 iNV 
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eo 53 ( Op2  13yi3 IOS. O22 
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APPENDIX D 


NUMERICAL INTEGRATION TI 59 PROGRAM 


A. INTRODUCTION 


In order to perform the numerical integration of on(a)da, 
a standard Texas Instrument's program was used [Ref. 25]. 
This program performs the integration by using Simpson's 


discrete approximation based on the following expression 


x 
n on 
—— J Siecle eedeg en oh ie Ok ates coed eee tt) 
O 


where 


f(x) must be known at n+l equally spaced points (f£ -f£,)- 
=o 
h = , X_>x 
" n 





O 


n+7 = number of data registers available 


n=number of subintervals = 2, 4, 6, .«..-. 


B. PROGRAM INSTRUCTIONS 
The program is taken from the master le br atayee proce ait 


package by using the code 2nd Pgm 10 on the calculator. 


a 








The user's instructions are as follows 


Procedure Enter Press Display 
a 
Enter data Subintervals A nt 
Enter data h B ht 
Enter data Function values 
0 Cc 
£ R/S ft 
Oo Oo 
1 
fy R/S fy 
i 
ae R/S fn 
Calculate = D Te 


i these values are printed automatically if the calculator 


is connected to PC-LOOA Print Cradle. 


The program steps are listed below giving location (LOC), 


key symbol (KEY), and comments. 


000 76 LBL 019 99 PRT 038 06 6 
O01 114A 020 92 RTN 039/545) 
@@2 53 ( O21 76. Let 040 42 STO 
G03 50 1xl 022 52 EE OA One Ot 
004 42 STO 023 00 0 O42 Ss 2a 
005 05 05 024 35a 043 98 ADV 
006 55 : 025 92 RTN 044 92 RIN 
987 O2 2 026 76 LBL 045 76 LBL 
008 54 ) 027° 7278 046 50 1xl 
009 42 STO O28) 42 7SR© 047 76 ST* 
010 02 02 029 03 03 048 O01 Ol 
O71 22 INV 030 99 PRT 049 32 xSt 
fe | OO ONT 031 92 RIN OSOmmO 1. 1 
Oils 29 cP 032 76 LBL O51 44 SUM 
014 22 INV Ogee ose 052) 01-01 
O85 67 EQ 034 53 { 053 32 xSt 
016 52 EE 035° 624-CE 054 99 PRT 
Oa7 «643 REL 036 85 + 055 92 RIN 
ome 05 05 Os ete 056 61 GTO 


Pao2 








1x1 
LBL 


LZ 


SUM 


RCL 
03 


PRD 
04 
REL 


ADV 
PRE 
RTN 


APPENDIX E 


TURN RATE LIMITATION TI 59 PROGRAMS 


fe) 6 LE RODUCTION 


The turn rate limitation calculations are divided into 


two parts with separate programs. 


Be. PART I PROGRAM 

For a given initial detection range and a given ship 
speed, Part I performs a sequential calculation of OarkKye 
One K,, 0% and K* for each of the maximum turn rates investi- 


gated, together with the "no limit" conditions (Vo) ur and 


Kor. based on 





On = on (BL elmax ; <99° 
Ss 
O 
o =—asQ = On 
Pp 
; : V,(1l+cos on B) 
BCU P-2 
ee" a ,B) 


o* = gos (P-2) 


2,1-P/2 


V, (1+P/2)° [1-(P/2) 


—- 1... 
? | max 
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POLE sp-2, the solution for * and K*¥ are not valid. 


Purchermore, for ie (Vo) no solutions are valid for any 


orf the quantities. 


The program is based on the fixed torpedo speed of 


< 
il 


rT 35 kts (18 m/s) and an initial detection range of 
T =e Ry = 250) i. 


pae wiser instructions for the program are as follows: 


Prodcedure Enter Press Display 
Enter data Maximum Turn Rate A lo|__ + 
max 
Enter data Ship speed B Vina 
Calculate 0* l 
and K* -- C o* ,K* 
Calculate $6 i 
A nee) o 
and ?p D Al B 
Calculate K E Ke, 
alculate —— ’ 
and K, - a 
al. 


These values are printed automatically if the calculator 
is connected to PC-100A Print Cradle. 


A listing of the program steps follows giving location 


(LOC), code (COD), key symbol (KEY), and comments. 


O@O 76 LBL Ol2Z 62 OF 024 54 ) 
001 114A 013 00 00 OZ>m 42 °STC 
002 42 STO O14 0s..3 (iZom Osos 
Hes O1 01 OUS, -0s2 OZ ymeae 2 Or 
ime 99 PRT OFG 62eOP 028 06 06 
Ug 91 R/S 017 04 04 GZ co Jerr 
006 76 LBL OUese OO sum 62 OP 
Gay 128 O12) Ole 031 00 00 
008 42 STO 020 08 8 OZ 02 2 
oom O02 O02 02a Sah = 033 06 6 
iso «699 PRT 022 43 RCL 034 03 3 
er 62> «CLR 023; 02802 OS Saeed. Al 
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510) 7) 


309 
Sa0 
S11 
Biz 


314 
a5 
SLE 


318 
Bug 
320 


Se 
325 
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B22 
S16 
Bay 
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Big 
350 
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333 
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242 $75 


Bas 02 2 
344 54 ) 
345 54 ) 
346 54 ) 
347 54 ) 
oeo 60682 «CSTO 
3aS 09 09 
eae 69 OP 
pa COG «(06 
fee 91 R/S 
bee O02 2 
eae 54 ) 
Bas 54 } 
e2>o 54 ) 
Say 864) 
ees 42 STO 
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C. PART II PROGRAM 
Given K from Part I, Part II program performs the probe 
calculations of the corresponding initial approach angles 


i 
(oo) ar (06) Br ($4) , and (O0) wx based on the relationship 


| eon 
E(o,) = Ka pie wi lSinlooda pe wr) 


e 


-ry[(1+(cosé.) a Be wy) }=0 


The program used a fixed torpedo speed of Vn = SS) kts le m7s) 
and an initial detection range of eee Ry = 250 m. The purpose 
of this general probe program is to locate roots of the given 
function y = £(o,) to evaluate the slope of the tangent line, 
and to find the maximum and minimum points on a graph. We 


will only use the first feature. 
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The program requires a subroutine for the function to 
be investigated. This subroutine starts at program location 
140 and 1s located at label 2nd C. The user's instructions 


for the program are as follows: 


Procedure Enter Press Display 

Enter data Starting Value 
ee om 2nd D ve 

Enter data The increment Ao, 2nd E Ad, 
Calculate 

Set Abe) == A £(o +46.) 
Calculate 

£(¢,-A9)) -- B £($,-46,) 
Display current 

value of ®% -- E %% 


If the value of Ad, 1s chosen too large, Ad, may be 
replaced by 4o,/10 by pressing label 2nd A. Similarly, if 
a larger value of Ad, 1s required, Ad, can be replaced by 
1049, by pressing label 2nd B. A listing of the program 
steps follows, giving location (LOC), code (COD), key symbol 
(KEY), and comments. 

Associated with the subroutine, it should be noted that 
the value of K,L, and p are entered separately in the 


memory locations 10, 11, and 12, respectively. 
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APPENDIX F 


FIGURES 
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Fig. 2. Location of Weather Stations Along the Norwegian 
Coast. 
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Fig. 6. Resonant Bubble Density at 38 kHz as a Function of 
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Fig. 12. Frequencies of Wind Directions in Percent for 
Stations Along the Coast from "Hillesoy" to "Ona." 
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Fig. 13. Frequencies of Wind Directions in Percent for 
Stations Along the Coast from "Myken" to "Furuholmen." 
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Fig. 24a. Geometry at the Sea Surface Scattering. 
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Fig. 24b. Specular Scattering Geometry. 
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Fig. 25. Rough Surface Shadowing Geometry. 


ILS) Fy 
















= source depth 
h,. = receiver depth 

Od = half beam width of receiver 
W = Wind speed 

Rp 


Varying: detection range 










Geometry : 
- 9 


From: © <*> 


















a mae! 
° Vv 
-R, -erf(v) 


© Crfe(v)#l-erfy 
° B 
- 5(Q) 


£$2>_ AS(O)R™ 
‘PP? R2gnd$'?> 


. Dy. 
ch tal IO leg Ia 


Te es 










rhpinect 
ore ae — Llpoirect 


Fig. 26. Calculation Scheme of the Surface Scattering Effect. 
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Fig. 28. Geometry of the Subsurface Propagation Model. 
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Fig. 30. Scattering Pattern. 
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Fig. 38. Pursuit Homing Geometry. 
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Fig. 40. Hit Criterion Geometry. 
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Fig. 41. AOB-limitation Versus Ship Speed with Parameter 
Maximum Turn Rate. 
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